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SYNOPSIS 

\ 

DESIGN OF A 120 MWt REACTOR FOR A SURFACE WARSHIP' 

by 

LT. ARUN J. RAO, IN 

Never before has the strategy of maritime warfare 
diversified so much as with the ingress of the nuclear reactor 
onto warships. Subsurface or surface, the nuclear propelled 
warships are far superior both tactically and in performance over 
their conventional counterparts. With the fast dwindling number 
of warships due to their esclating costs, the possession of ! 

nuclear-powered warships for maintaining ’the same effectiveness 
and punch is becoming increasingly inescapable. 

Reactor design is a multi-dimensional complex problem^ 
embracing many fields of engineering. The warship reactor design : 
is more intricate than the land-based ones at it involves 
certain additional specialized problems like ship-motion effects, 
manoeuvring, vibrations, etc. The physical parameters of the j 
core is restricted mostly by thermal considerations. Hence j 

f 

the Thermal hydraulic design requires careful attention to 

I 

prevent core melt-down or coolant bulk-boiling. The Reactor 
physics design involves evolution of core enrichment and certain 
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other parameters to result in optimal design. Knowledge of 
the neutron flux distribution is also required. Here, two-group 
theory for heterogeneous core lattices has been used. The 
Shielding design involves meticulous calculations to produce 
the smallest possible reactor shield, both in size and weight 
which is of paramount concern in a warship design. 

Control of a marine reactor is different from 
the land-based ones. Due to stringent manoeuvring requirements. 
Power excursions from 20% to full-load require automatic 
control with a high power increment/decrement rate of 3 °/^/ . 

The system has to be prompt acting, stable at all times 
of operation with no unallowable overshoots. Finally, when 
the reactor goes on the ship, certain other problems like 
ship-motion effects, vibrations, etc., come in. 

Due to the immense calculations involved, any 
modern reactor design process is impossible without the aid 
of a digital computer. Here, 3 computer CODES have been 
developed, one each for thermal hydraulic. Reactor physics 
and Shielding designs. The CODES & their results have been 
presented separately. 



The different design processes discussea till 
now are inherently interactive among themselves and cannot 
be considered isolatedly. This nas been highlighted in the early 
part of the thesis. It has, however, been the earnest t-nd.e^v^ur 
of the autnor to ensure an overall non-conflicting design. 
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INTRODUCTION 


In the modern era of fast depleting oil,, coal 
and other forms of fossil fuels, nuclear energy takes a 
leading seat as the future global energy source. The. nuclear 
reacfor is virtually an interminable source of energy. Tapping 
of this energy is an engineering feat by itself. Though 
marginally more expensive,, it is only a matter of time in this 
energy-crisised. world, before the reactor, establishes itself 
as an inevitable source of energy. All this, and it is a mere 
42 years old 


Today the cost of a nuclear reactor is 
prohibitive. Why then nuclear propulsion? The answer is 
straightforward - No repeated- refuelling for the ship;, which 
opens an, entirely new era at sea.' Where economic consideration 
is counterbalanced by military requirement,, nuclear propulsion, 
for -Naval, strips become essential. For the merchant ships, 
however, this Is not so, as money is the key factor for 

trad?. Consequently, we have ojjjy a handful of Nuclear merchant 

whereas the number- of . ,, , 

propulsion) exceed three hundred. P world ^ith nucha? 

SCOPE OF THE - THESIS 


A compact low weight system is of ..paramount 
importance in plants designed for. ship propulsion. Marine 
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reactor design by itself/ is a very vast topic merging land- 
based reactor design with certain specialized aspects like , 
marine environment/ light weight/ limited space/ etc* Our 
goal is to design a reactor system for a frigate/destroyer 
class warship* A drawing board design involves a vast amount 
of data. } ti|ne and generation of numerous highly essential 
computer codes. This is definitely not being attempted. In 
the following pages# certain basic aspects of marine reactors 
have been discussed and a skeletal design attempted to 
highlight the varied & complex problem of marine reactor 
design. It is beyond the scope of this thesis to cover all 
topics envolving the design# be the skeletal stage itself.. 

All important aspects have#, however# been covered. 

LITERATURE REVIEW 

Most of the information regarding Naval reactors 
ic classified. Though information on reactor sof . merchant 
ships are easily available# they ..don't satisfy the stringent 
measures of manoeuvrability , .high power- to- weight ratio. . - 
and the limited space available in a warship where economic 
considerations fall into the background* Yet# a fair idea 
on problems encountered in marine reactor design#, is available 
from them.. International symposiums on Nuclear ships have , , 
been held in (1960)# Hamburg (1971) .and 

Madrid (1972). 
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A subcommittee on military applications of the 

U.s. has presented the tactical and effective use of Nuclear 

Propulsion for Naval Warships It argues that, though 

economically dear, when faced with the threat of submarines 

armed with either torpedoes or cruise missiles* or bombers 

or surface ships armed with cruise missies, the all-nuclear 

task force (consisting of an aircraft carrier and a few 

destroyers/frigates, all nuclear— powered) is less vulnerable 

to attack than a task force of conventional ships. The 

committee concludes, nuclear frigates for independent missions 

are far superior in performance than their conventional 

counterparts. John R Lamar sh^^ states that submarine reactors, 

?.35 ^ 

being very small and compact, use fuel enriched in u to .. 
over 90 w% A.M. Petrosy.'ants >©^ 4 .cidating problems of 
nuclear science Si technology, discusses Russian, ships with 

nuclear engines - 'Lenin' & 'Arktika' (both ice-breakers) j 

I 

& their vastly superior performances over diesel icebreakers. ! 

i 

f 

Extensive study & analysis of the world’s first j 

(5) . • ^ 

nuclear-prapelled merchant ship 'Savannah , in, its sea | 

environment, has thrown light over the problems faced by, 

Marine reactors. Analysis of its core performance, has. been, done 

by WJ. Gallagher and Kaiz ^ . The Japanese Nuclear ship 'Mutsu' „ 

( 6 ) 

has also been discussed in detail by S. Sasaki Off-late, 

L. Chinaglia^ 2 * of Fiat, Italy has presented a new design of 

, 

f 


[■ 
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a PWR ,for ship propulsion. The inherent design interactions 
among various design aspects like thermal hydraulic design, 
nuclear design, shielding design, etc. has been well 
illustrated;-, by ,RA Knief ( 13 * . 

El-Vakil ( 17 J deals with PWR thermal hydraulic 

( 26 ) 

design. Weisman ' describes at length the various hot-spot. 

factors concerning the thermal design. Regarding .reactor 

( 22 ) 

physics analysis, Meem , with his 2-gp theory, gives 

sample calculations for .reactor designers, that has been 

( 21 ) 

extensively referred. The Nuclear Engineering handbook - . 

^by Etherington has been used for related figures & 

references as and wl^en mentioned. IAEA reports Vienna, 

(27 ) 

1964 , highlight the reactor shielding design .and problem. 

B. Chinaglia^®^ of Fiat, Italy, has also worked on shielding 

studies & control of ship-based gaseous activity release. 

(29) 

J. Brown & A. Jackson . have carried out a study on 
shielding of Marine reactor installations. S.J. Ihunell^^ 
has presented a comparative study on neutron fluxes in a I 

metal-watep shield in merchant ship reactors. 

[ 

A Marine Reactor control & in-core instrumentation f 

I? 

system design has been done by F. Basile, R. Merino & Q, Selviaaa. 1 

l 

f 

wherein they have pointed out the differences between land* 

based reactor, & marine reactor control. The control scheme j 

( 2 ) ! 

used in the thesis is one evolved by Chinaglia . of Italy i 
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A detailed treatment of controllers has been illustrated 
( 33 ) 

by £4A Schulte-: * this has been used, however, 

after a study of non-linear control system. References for 
other aspects on reactor control have -been Lewins^ 3 ^ and* 
Gladstone & Sesonpke^ 12 ^ and D. L. Hetrick^ 36 ^. 

Specialized marine problems like ship-motion 
effects on steam-generator performance, collision barrier. 

6c grounding protection barrier for- ship s^re actors, and 
safety of crew have been studied, respectively, by V. Jecob 
and G. Selvegg^ 3 ^^, W. Jager & H. Lettnin ^ 33 ^ and L. Chinaglia 
& P. Rucci^ 32 ^ . 

1.3' THESIS ABSTRACT 

The thesis starts with a brief history of marine 
reactors. The tactical aspect of nuclear— powered vessels 
over their conventional counterparts is discussed- this may * 
m^ke the difference between winning or losing a battle which ... 
is vital for any fighting force. The salient features of . 
ship-board reactors over the land-based ones.g^ discussed. 

The modern trends of marine reactors aj; . e a Iso elucidated. . . 

Chapter .3 highlights the. design process. The . 
different aspects of reactor design like thermal hydraulic 
design, nuclear • design, shielding design, etc., are inherently, 
interactive among themselves and hence cannot be independently 
considered; this is highlighted. A reactor core configuration^, 



7 


is then evolved and the systems used mentioned. 

Core thermal hydraulic design is developed in 
chapter 4. A linear power rating of 11.0 KW/m is used. 

Since the power generated by a reactor is essentially resti -Ctg^ 1 
by the effective heat removal, this portion of the design 
is carefully considered as we desire the most compact and 
light-weight reactor for ship-board installation. A computer 
CODE, so vital- .for any thermal analysis, has been developed 
and used to arrive at optimum reactor Parameters. With this 
knowledge, we next move over to the reactor physics design 
(chapter 5). Two-group theory is used and the core studied 
for criticality using 2 regions (water as reflector). A 
computer CODE has been evolved which computes the effective 
multiplication constant (K^) , the fast and thermal fluac.es 
in both the regions, as well as the various 'poison 1 (unwanted 
neutron absorbants) activities. in the core. The amount of 
control poison required for reactor control is also evolved. 

The control design hence would chiefly depend on the parameters 
obtained in this chapter. Care has been taken in the COPE 
that temperature and,, pressure dependencies of all parameters 
ar eopscwau&ider ed . 

We are aware that the overall reactor weight 
and size has to be the minimum possible for a marine reactor. 



8 


Bui it of this weight,, however, is taken up by the biological 
shields-. Consequently,, re.actor shielding assumes considerable 
importance and is analyzed at length in chapter-6. Biological 
shielding is mainly required for the fast neutron and gamma 
fluxes emerging from the shield materials due, to thermal 
neutron capture. The biological dose rate from the secondary 
shielding is to be restricted to 0,75 mrem/hr. This has been 
achieved after the formulation of a shielding CODE. In 
addition an approximate weight calculation of the reactor 
(with shield materials) is demonstrated. 

The, control of a marine reactor is greatly 
different from a land- based one. This is because of high 
manoeuvring demands and a requirement of controlling the 
reactor automatically . from full- load to part-load (20%),,, 
Chapter-7 elaborates these aspects and evolves a reactor control 
system. Careful stability analysis is also done from full-load 
to part- load using a two temperature feedback for the 
reactor. Certain other aspects of stability are also attempted 
after the calculation of the temperature coefficients of 
reactivity as applicable. A reactivity balance of the core 
using the CODE of chapter 5 is displayed. The control policy,.- • 
used is elucidated subsequent ly. 
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The integrated aspect of the reactor with 
the ship increases the design complexity further. Its effects 
are of practical significance to the operator on-board. 

Chapter-8 embellishes this integrated view. The ships threc 
dimensional motion at sea has its effects on the steam 
generator performance. Collision barrier and grounding protection 
structure is required for the reactor compartment. 

Chapter-9 concludes the present work 
illustrating the design that is ultimately arrived at. Other 
aspects of marine reactor design beyond the scope of this 
thesis are mentioned and suggestions for further study 
given. All references made in the current task are listed 
in the bibliography. 


* * * * * 
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In 1942, Fermi & his co-workers made the world *s 
first reactor critical at the University of Chicago, However, 
not until 1957 was a propulsion unit designed to operate on 
nuclear fuel. It was a small highly-enriched Pressurized 
Water Reactor (PWR), installed on board the Submarine USS 
NAUTIHJS. As of today, there are over 300 ships propelled 
by nuclear energy. Most of these, however, are warships 
where economy is not the deciding factor in the choice of 
propulsion system, 

2 . 1 BRIEF HISTORY 

2.1,1 NUCLEAR MERCHANT SHIPS 

The history of nuclear merchant vessels have never 
been bright. Nuclear power in merchant ships was used for 
the first time on a 16,000 ton Russian ice-breaker 
in 1960. The decision was based not. on economy but efficiency 
over the prevailing diesel electric system. In Arctic areas 
of heavy ice, nuclear ice-breakers with a 50% higher 

specific capacity*, can operate most of the year; but the 
conventional ones can make way only for 3-summer months 
due to danger of running out of fuel. The second ship was 
‘SAVANNAH' , a U.S. Cargo-liner commissioned in dec *62 
(no longer in active service). It's purpose was to test 
safety of nuclear plant aboardship & entrance of nuclear 

*Specific capacity - power-to-displacement ratio. 
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vessels to civilian ports. 

The 'OTTOHAHN* of ERG is a small scale experimental 
ship not planned for profitable operation (1968 vintage). 

Her purpose is to furnish all types of information theorio .cal, 
practical, economic & scientific with regards to nuclear 
propelled merchant ships at sea & refute arguments about their 
safety. The Japanese, on similar lines as ERG, launched a 
10,400 ton specialized nuclear cargo ship in 1970 (MUTSU). 
Russia, finding the need for a second nuclear ice-breaker 
commissioned 1 ARKTIKA 1 (displacement 23,460 tons) in 1975. 

The reason for this bleak growth of nuclear ships 

arc plain 

(i) Cost economy - Capital cost of nuclear ships are 60% 

higher than conventional. The fuel cost 
is also higher. Till such time fossil fuel 
becomes prohibitive enough to counter- 
balance the above two reasons, nuclear 
ships will not be easily accepted. 

(il) Uncertainties-of investment cost; the future development 

in this area is hardly predictable. 

(iii)Safety problems- Nuclear ships are not easily accepted 

by most of the civilian ports; psychological 
fear of public is mainly responsible for 


this 
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2. l t 2 NUCLEAR WARSHIPS 

, A 

, , Nuclear Energy toox birth In warships in the 

form of a submarine reactor for 1 NAUTILUS* , a U.S» Submarine, on 
17 Jan '55. She generated 15,000 SHP, having 2 shafts and 
displacing 4040 tons. Initially- she steamed submerged -from' -Key 
West, Florida, to New London, a distance of 1397 sea miles£ iis 
She was refuelled after 2 years and 62,562 miles of steaming. 

Her success gave rise to a 'gold rush' and the US Navy ordered 
a pack of. submarine reactors from .Westinghouse Corporation. In 
i960, Russia started her submarine fleet. 

•ENTERPRISE', a 89,60C » tons US Attack carrier., 
was the first surface ship to have nuclear propulsion ,Cin 1961,),* • 
She has 8 PWRs on 4 shafts, giving her a power of 28Q # 00Q ..SHP. 
and a maximum .speed of 35 knots**. She has had her core changes! , , - ; 
thrice, with the present feed giving her a sustainence of 10 
to 13 years. In. 1962, ' BAINBRIDGE ' , a U.S. Nuclear-Powered 
Guided Missile. cruiser was the first cruiser to be nucle.ar- 


propelled (displacement 8580 tons). 

Changes in. ship reactor -technology have 

also .been*. ..incorporated to make -ships • more versatile. end~: : — - 

independent... The newer -US. aircraft carrier^. ^NIMXTZ', (l^pyached,_- 


in 1.9,75) has only 2. reactors contrary to her 



’ENTERPRISE* witir-S; this- generates a ..power of ^260‘,000 ’SHP ©r 


4 shafts,., giving her ' 93., 400; _ tons a maximum speed, p £ . .35 - knot 


She is the largest warship ©v@r. . . 


* I se,e mile (or nautical mile) = 1.83 kms (2000 yar^a) „ 

** 1 khow = 1 nautical mile/hour. 

\ - ■ 
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built. She can steam an estimated 800,000 to 1 million miles 
before refuelling. Interestingly, the energy equivalent of 
her cores is 462 mega gallons ( i.<&0 million tons) of fuel 
oil- this is where nuclear energy scores over the conventional 
form. 

The largest submarine ever built, however, is 
a Russian one, a 'typhoon* class Ballastic Missile submarine 
(SSBN), launched in 1980. She displaces 30,000 tons dived 
with speed exceeding 24 knots and having 2 nuclear reactors, 
one per shaft, supplying a total of 120,000 SHP; she carries 
20 SLBMs ( Submarine Launched Ballistic Missiles). Presently,, 
nuclear submarines baas outnumber nuclear surface ships. 

However, many surface ships both cruisers and destroyers 
are increasingly becoming nuclear-propelled , mainly due to 
tactical reasons. 

Due to the prohibitive cost & superior technology 
for marine reactors, nuclear warships are expensive toys 
handled only by a few economically advanced countries. However 
with the Indian Ocean swarming with subsurface &>_ surface, 
nuclear & conventional warships, possession of nuclear warships 
to effectively safeguard our exclusive economic zone . & 
sea-fr ontiers, is becoming increasingly inescapable. 
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2 . 1,3 FUTURE TRENDS 

For merchant ships greater than 10,000 tons 

, (23) 

displacement, preliminary studies indicate nuclear Power 

to sc .re over conventional power when operated for long 

durations. This might be used for constructing huge nuclear 

container ships in future when nuclear energy gets accepted 

in the merchant marine field. The smallest nuclear ship ever 

built is a US nuclear -powered ocean engineering & research 

vehicle 1 NR-1* (1969), a submersible, built by USAEC*, 

Knolls' Atomic Power laboratory. The vessel displaces only 

i i • 

400 tons, having dimensions of 136.4 x 12.4 x 14.6 with 
2 propellors & carrying 5 crew members and 2 scientists. 

Such vessels are very useful for conducting detailed 
undisturbed experimentation on the continental shelf over 
long durations. 

Now we hear the concept of NUFIS - Nuclear 

Floating Island, brought forward by Japanese scientists. 

Because of small area & dense population, the siting of 

Nuclear Power Plants in Japan is limited. Hence, the feasibility 
(9 ) 

of NUFIS has been studied and the project is about to 
commend^. It is to be located 20 kms offshore at a sea-depth 
of 150 m. It is to have a 1200 MM^ FWR Plant on a 140 m x 140 m 
floating platform with a displacement of 298,000 tons. The 
power plant &. the associated facilities are to be situated 


4 
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on it. Hence, the marine reactor technology will have to be 
applied to it. The technological feasibility of a floating 
centre for fuel fabrication, reprocessing, etc., is also 
being thought' - :. 

2.2 TACTICAL ASPECTS OF NUCLEAR WARSHIPS 

It was stated that the elimination of repeated 
refuelling opens an entirely new era at sea. Let us explore 
the tactical boost this gives. 

Nuclear propulsion enables warships to traverse 
across seas at high speeds without parsing to refuel and 
without replenishment. With the escalating cost of warships, 
many many times that of equivalent merchant ships, the 
number of warships in the fleet is bound to dwindle. Hence, 
'nuclear propulsion becomes prominent to maintain the same 
effectiveness with lesser warships, with the fast depleting 
oil sources^ in future, nuclear propulsion is going to be<pome 
a necessity. It gives a n additional tactical f lexibility 1 and 
opens new courses of acti-.n to the operators of the Naval'., 
Forces which would otherwise be unavailable. A nuclear 
powered naval force will be able to achieve many more \ 
productive results in any operation; such a force would) be, 
in addition, less vulnerable to the action of the enemy, U 

*' v 

/ 

\ 

\ 
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2.2.1 DRAWBACKS OF CONVENTIONAL SHIPS 

Here, mobility is limited by the need for refuelling 
at sea. Enemy action, bad weather & other uncontrollable 
' forces directly impact upon refuelling ability which is not 
only time consuming but also requires the services of a 
considerable © 43 . tanker fleet/train. The Fuel oil carried by 
ships being inflammable, in addition, poses a fire .hazard 
onboard. A nuclear fleet, contrarily, is freed from the 
constraints of tankers & base support; hence it has increased 
speed, reliability, range, less fuel requirements & more 
payload . 

2.2.2 NUCLEAR SUBMARINES 

In the modern warfield, no other weapon has proved 

to be such an imminent threat as the nuclear submarines. As 

we see from Table 2.1, the super powers have maintained a 

parity on the number of nuclear submarines. Admiral Moorer 

( 4 ) 

of US Navy has stated that he would rather have one nuclear 
powered . tactical cruise missile submarine (SSbn) than 5 of 
its conventional counterparts. Submarines are a particular 
threat because of their potential for approaching to within 
weapon ranges without being detected. 

When their batteries run down, diesel submarines 
have to periodically ‘snort? 1 . The quantity of fuel they carry 

** Snort* - come up to near sea— surface level to to run diesel 
engines for charging batteries that are used for 
propulsion underwater. 



limit their edurance. The nuclear submarines, ostensibly, 1 
never require to surface and can continuously move at high 

* 4- 

speeds of 20 — 30 knots ( almost double of diesel submarines! )• 

The modern nuclear submarine being more quieter*, have a 
high probability of going undetected by enemy sonars. Hence,'' '' 
they are much more effective in penetrating anti-submarine 
defences than diesel submarines. 

2.2.3 SURFACE SHIPS - THE ALL NUCLEAR TASK FORCE 

A naval task force consists of an aircraft 

carrier with her modern air-wing, and her escorts of 3 to 4 

destroyers/ frigates to provide additional protection against 

surface, sub-surface and air threats. The capability and 

punch of the force is restricted by the endurance of each 

of its ships. Hence, for maximum effectiveness, all ships 

in it have to be nucl ear-propelled . When faced with the 

threat of submarines either armed with torpedoes or cruise 
missiles, or bombers, or surface ships armed with cruise 
miaggi&ea, the all-nuclear task force is less vulnerable to 

attack than a task force of conventional ships which must 

show down to refuel in an extended encounter. It is therefore 

much more difficult for the submarine to get in a favourable 

position to attack the all-nuclear force. Further, the 

logistic support forces upon which the conventional striking 

forces are more dependent, are even more vulnerable to attack the 

* , quieter* - Vibrations of hull due to inner machinery is 

of lesser magnitude, making passive underwater 
detection more difficult. 



warship^. 


The tactical advantages of nuclear propulsion for 
ships of a nuclear task force can hence be summed up as 
follows 

They can move where required at top speeds without 

establishing logistic support before departure. 

Their virtually unlimited range enables them to use 

their superior speed to arrive at crisised situation 
much earlier than their conventional counterparts. 

When they do get there in a less vulnerable condition/ 

they have far greater staying power, not only in terms 
of their ability to stick on without fuel, but also, 
significantly, in terms of aviation fuel & aviation 
ordanance for their airwing. For example, ^dbs 
' NIMITZ 1 has 50% more aviation ordinance & twice as 
much jet fuel than equivalent conventional carriers, 

It can continue combat operation without support 

until aircraft fuel, ammunition, food or other 
supplies are exhausted. 

It can disengage if required to retire at top speed 

to less vulnerable sources of aircraft fuel & ammunition 
replenishment. 



2.2.4 NUCLEAR PROPULSION FOR FRIGATES 


By virtue of posses ing their capability for high 
sustained speeds and uninhibited endurance, nuclear frigates 
are uniquely suited for assignment to independent tasks with 
an effectiveness no non-nuclear ship can match. This can mean 
difference between winning or losing engagements. In independent 
operations, it matches the endurance of any nuclear submarine. 
Since the aim of the thesis is design of a reactor for frigate / 
destroyer, let us explore its tactical aspects a bit deeper. 

Nuclear frigates combine the essentially unlimited 
endurance of nuclear propulsion with the capabilities of the 
most modern anti-air and anti-submarine weapon systems. Their 
nuclear propulsion results in major improvements in offensive 
& defensive capabilities 

greater attack effectiveness due to increased mobility 
& freedom of independent action; ability to transit at 
high-speeds and thus arrive in attack position earlier 
& ready to conduct operations immediately without pausing 
to refuel; ability to be on attack station a higher 
percentage of time; the freedom to extend attack to 
larger areas, 

reduce vulnerability due to freedom from dependence upon 

replenishment in contest areas; and ability to avoid 
submarine attack by transiting at high speed Sc . using 
evagive transit routes. 

significantly reduced need for sea-based & land-based 

logistic support forces. 
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,-incr eased propulsion plant reliability duetto the 

higljer standards required in the design con struct ion., of 

nuclear propulsion plants. 

increased overall effectiveness of the entire fleet 

by reducing the total number- .of comhat ant. ships,. required 
to protect , oilers . 

2.2,. 5 WORLD 1 S NUCLEAR FLEET STRENGTH 

Table 2.1 lists the nuclear strength of the world 
navies as of today. Though nuclear submarines for outnumber 
nuclear surface ships., we can see that,, the numbers of the 
latter has more than tripled in the last. 10 years and many 
are being built* 
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.-MARINE R5 ACTORS 


LAND BiwSED REACTORS 


The crew (& passengers) 
constitute a 'captive' audience; 
they have no place of shelter 

in case of accidents. 

In such an event, - 

people could be evacua- 
ted to safer areas. 

Most likely place of accident 

is in the least favourable 

location - the port areas. 

This problem is 

non-existent , . 

(ii) DESIGN 


There are many more start-up 

and shutdowns as well as widely 

varying needs for power in man- 

ctauvering. This results on a 

greater stress & strain on the 

nuclear machinery (to be taken 
care of in .design. ) 

The design need not 

be that stringent. 

i 

Due to space limitations, the 
machinery of working piantss is 
much more crowded than in 

stationary plants. 

No such limitations.. 

. Due to weight & space limitat- 
ions, the shops, spare-parts 

and maintenance facilities are 
• 

kept to a minimum. 

- do - 

Design has to cater for specialist 
aspects on ships like vibrations 
(of propeller & machinery) , sea 
conditions (roll, pitch, yaw) and 
limited space and weight. 

Earthquake protection 
is required for reactors. 
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lA^-BASSD REACTORS 

11 . 

Irrespective of economic : 

Materials chosen should 


considerations/ only the best . : 

be such as to minimize 


materials are to be chosen such j 

capital & maintenance 


as to achieve highest possible 

| 

cost while still being 


specific capacity. 

effective. 

12. 

Reactor control should be such 

aa to make it inherently safe* 

in the event of a tsunami or 

sinking of the ship. ! 

1 

, No such restriction 

13. 

j 

Propulsion system design should 
cater for an emergency conventi- 
onally powered propulsion system; 

No such condition. 


coinpat able with the main machinery/ 


in the event of a complete break- 
down. j 

1 


14. . 

! 

In warships/ adequate shock j 

No such condition 


protection for reactors is required required to be 


against firing of onboard gun® 
and missiles? also required .is 

(i) protection against shelling 

(ii) grounding & collision 
protection against external 

accidents. 

catered for. 

15^ 

Only PWR ' s are used as marine 

All types of reactors 


reactors. 

are possible on land. . 
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2.3,2 BASIC . ggaifiN CRITEPIA FOR -NUgl&AK K&SSHIRS. 

Threats to crew health & safety do not arise 
out of normal operation of nuclear ships but rather du,e 
to the possibility of internally and externally caused 
accidents to them. To minimize these, a variety of special 
safeguards in design, construction & operation are required. , 
These will be briefly elucidated. There are also, a variety 
of recommendations applicable to nuclear ships; these, 
however > won 't be discussed here but can be found in 
references (7) & (10). 

Considerations such as the following lead to, 
certain basic objectives that must be achieved in the design 
of a nuclear warship; 

- There must be continual knowledge of the level of radio- 
activity in, all normal effluent streams,,. & of the levels . . 
at various locations around the ship. 

- There mu.st be complete control, of .all normal effluent 

st r e ams , ,wit h storage capacity sufficient to permit detailed 
release of effluents under unfavourable circumstances 

- There -must -be.- an extraordinary degree of conservatism/ 
redundancy & dependability in the design of the reactor, 

so that the likelihood of .accidents will be. reduced, -tp the . 
Ipwent practicable level. Following -aspects are important s - 
Marine environment like roll, pitch, yaw< 
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<*/ PrQpeller & machinery vibrations of the reactor 
dompartrnent . 

^ Proper reactor containment design for safety, 

^ Collision & grounding barrier protection for 
reactor compartment. 

/ 

Protection of reactor against shocks impressed 
by weapons firing onboard. 

jJf 

The ship should have highest possible specific 
capacity. 

- There must be devices & safeguards to minimize the hazardous 
consequences in case serious accidents should occur. 

- There must be an adequate conventional emergency propulsion 
system in case of complete main-engine breakdown. 

Thus the importance of nuclear propulsion for 
warships as a means towards strategic mobility is obvious 
to all who pause to consider the matter. It will remain as , , 
until a superior form of ocean propulsion for the weapon 
systems is produced. 
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It is well known that only PWRs are used as mobile 
reactors. In this chapter, this has been justified. In 
addition, the design methodology, the fuel assembly design 
and the . steam-cycle used has also been elaborated. 

3.1 CH OICE OF FEACTOR 

"S&AWGLF" , * the second nuclear -powered U.S. Submarine 

was fitted w.ith a liquid metal (JNlak) reactor called the 
Submarine Intermediate Reactor (SIR) in 1956. During dockside 
trials, steam-leaks deveoped due to Nak alloy coolant entering 
the superheater steam-piping. After repairs' she was put to sea. 
However, due to recurrence of the above leaks in Dec ’58 , 
the reactor had to be replaced with a PWR. 

The above experience was an eye-opener for the choice 
of mobile reactors. Alternately, if BWPJ were. used, they would 
give inherent system instability due to the reactor water 
level constantly varying with the ship-rmotion effects; in 
addition, for direct cycle the turbines require shielding, 
hiking up the gi$iss. .machinery weight which is a very undesirable 
phenomenon. 

Here, it is not possible to compare the various 
types of reactors due to space limitations. However, the 
salient features of a PWR over other reactors have been listed 
below:- 
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Advantages 


, - Moderator and coolant same. 

- Water easily available ana cheap* 

~ Water poses no serious long-term radioactive hazard. 

- A primary heat exchanger leak is not hazardous to 
personnel as. primary loop is shielded. 

- Reactor has a. high negative temperature coefficient 
making it safe & stable to operate. 

_ For moderate enrichment and closely— packed lattices, 
burnup is the highest. 

- Being pressurised, the reactor is compact and 
size a minimum due to high power densities. 

- PWR technology is very well known as it is. the 
most common type of naval reactor today* 

Pis adv.ant ages 

- Hicfti primary loop pressure required due to the low 
bpiling point of water; so pressure vessel should 
be very strong. 

-Poor thermodynamic ef f iciency-20-25% for mobile 
reactors. 

- As. ordinary (light) water has a large neutron 
absorption cross-section neutron economy is poor. 

- Reactor should be shut-down for reloading; 

( i.e. reloading not possible at sea). 
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i Despite its moderate disadvantages, the power-to- 
weight ratio, of a PWR is the highest; this is a highly 
desirable, .factor where compactness and light-weight is ,a chief 
requirement. All marine reactors operating today are hence 
only PWFs . 

3.2 THE DESIGN SCHME 

3.2.1 THE PROCESS 

This may be divided into the following categories:- 

(i) Core thermal hydraulics design 

(ii) Core Nuclear design 

(iii) Shielding Design 

(iv) Reactor control and safety 

(v) Materials Design 

Since each of the categories is complex by itself, 
it 4>s „ being treated separately in chapters that follow. 

The last topic ( v ) , however, has not been included separately, 
but sketchily covered throughout. Being a reactor for defence 
purpose, .no study on cost -economics has been done. . 

3. 2, 2. INTERACTIONS 

The integrated reactor design requires substantial 
interaction among often conflicting design principles of 
(i) to (iv) above, safety considerations like biological dose , 



DESIGN SPECIFICATIONS 
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liiT^-tation , prevention of core melt-down/ etc ? provide very 
stringent constraints. 

Fig 3.1 gives a pictorial view of these interactions 
in the present design. Nuclear, thermal-hydraulic and 
material design components result in composition, geometry 
and operational specifications that characterize the marine 
reactor system. Nuclear design begins with preliminary 
cross-sections data. For the calculatianal modules, flux- 
averaged cross-sections with appropriate modifications 
to the working conditions have been used. The flux, power 
distribution and burnup contribute to a final specification 
of enrichment. The thermal hydraulic design interacts with 
power distribution to develop a geometry and its related 
temperature and flew distributions, limited by material 
specifications. The nuclear design also yields the reactivity 
requirements required for the control design. 

The control system design for marine reactors is 

intricate. Whereas a power variation of only 10% allowed, 
in land-based P'WR 1 sj/arlation from 100% down to ' 

20% is required to be automatically controlled that too in 
much shorter times for marine reactors. This calls for a 
non-convent ion al control system design with careful stability 
analysis. This is dealt with in chapter-7. 

3.3 FUEL ASSEMBLY DESIGN x The first step in core design 
is to evolve the type of fuel assembly required 




Fig. 3. 2 

The fuel assembly 
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The basic features of the assembly selected is shown in 
Fig 3.2 (Ref. (7)). A 17 x. 17 square pitch layout has been 
chosen with 282 cylinderical fuel rods per assembly. This 
conf iguration represents recent designs that enhance linear 
heat-rate limits ( Ref. (13)). The taxonomy of the fuel 
assemblies is as shown in Fig, 3.3 . This symmetrical 
arrangement employs an open lattice that permits some flow 
mixing between adjucent units which is a desirable phenomenon.. 
The diameter of the core baffle is 7.071 times the width 
of a fuel assembly. 

A constant fuel enrichment has been chosen 
throughout the core. For a linear power rating chosen as 
11.00 KW/m (see chapter 4) and 120 MWt power output, the 
number of fuel assemblies required are 29 (8178 rods). This 
gives a core height of 133.39 cms. The core-radius is 
dependent on the fuel rod-pitch (which will be determined 
later) , 

6 1 

There are/control-rodf poison fingers per assembly- 

each finger being of the same length a$ the cose. It cpntains 
B^C ( Boron carbide) enclosed in stainless steel tubes. Their . 
concentration and size are to be determined. The reactivity 
control is provided by 29 Rod Cluster (Soptrol (RCC) assemblies 
(1 rod-cluster per assembly). They have been constituted in 
3 banks A, B and C as shown in Fig. 3.3. No burnable poison 
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or soluble poison scheme has been included. Each assembly/ 
in addition, has been provided with a central instrument thimble 
space meant for a movable in core flux measurement detector 

3.4 STEAM CYCLE 

The steam cycle used is illustrated in Fig. 3. 4. 

The 2-loop system is in vogue in .almost all working PWR's. 

Th=; primary system contains water <at 2000 psi (14 MPa ). 

Coolant pump regulates flow in the loop and pressure vessel. 

The presurizer maintains the system pressure within a specified 
range. ^,,2 primary coolant loops supplying to 2 u -tube steam 
generators of 60 MWt capacity each are chosen. 

■ The secondary loop generates steam at 500 psi 

(3.5 MPa) at 100% load at a quality of 0.9975. Regenerative. 

been 

cycles for the turbines have^/used to improve efficiency* 

Further details about the steam cycle have not. been attempted. 

3.5- INFERENCE 

We have now determined few of the reactor core 
parameters and the system we have to work on in the. succeeding 
chapters. We also know about system interactions. While . 
considering the various design processes in the subsequent 
chapters, these interactions have to be borne .in mind for 
an overall non-conflicting design. This has been the earnest 
endeavour of the author. 


*********** 
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4.1 INTRODUCTION 

A nuclear reactor is capable of operating at very- 
high power densities. However, the core must be operated at such 
a power level that with the best available heat removal systems, 
the cladding and fuel temperatures anywhere in the core do not 
exceed safe limits. A fuel element damage, otherwise, will result 
either in the transfusion of the highly radio-active fission- 
products into the coolant, or in core- fuel meltdown. Both these 
effects are highly uriJ.es ir able '■*. Hence, the core design is 
dictated more by the thermal limitations rather than any nuclear 
cons ider at ion . 

GOVERNING EQUATION 

The rate of heat removal, (Q) is given by the simple 
relation (p:369, (17)) 

Q = h* A *At- 

ov^al^^^Area of Tempr.diff. between 

ht. transfer heat transfer fuel & coolant, 
coefficient* 

‘h 1 is governed by conductive heat transfer through fuel elements 
and cladding, and convict ive- heat transfer through coolant. As 
the fuel is U0 2 , its thermal conductivity is fixed. So fuel-; 
cladding should have least resistance to heat conduction, in 
addition to low neutron absorption cross-section, workability. 
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ability to withstand high temperatures and withstand the corrosive 
action of the coolant fluid; some of these considerations are 
conflicting; hence a compromise between high thermal conductivity 
and high temperature operation necessitated the use of stainless 
steel (£d-347) as the cladding material. Resistance to heat 
conduction is reduced by keeping the cladding thickness to the 
minimum. Heat transfer by convection to water depends on the 
properties of water at the operating conditions# and coolant 
velocity. 

"A" is maximized by having as many fuel rods 
as possible. Since a large quantity of heat is generated in the 
compact core# the maximum possible heat transfer area should 
be squeezed into it. 

"/it"# finally# should be high; high rates of heat 
release requires low coolant temperatures; this conflicts with 
basic thermodynamic considerations as a good thermal efficiency 
requires high coolant temperatures. Hence# to, kepp coolant 
temperatures as high as possible and yet retain a high At in the 
core, fuel elements have to operate at the highest possible 
temperatures which are limited only by metallurgical & burn-up 
considerations. 

Hence# it is seen that reactor core design necessital 


a compromise between often conflicting nuclear# metallurgical# 
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structural , thermodynamic and heat-transfer requirements. To 
ensure that the core is safe under all operating conditions 
and transients in our main motto. 

4.2 DESIGN CHARACTERISTICS 
4.2.1 SALIENT FEATURES 

The main problems encountered in thermal hydraulic 

design are :- 

(i) assure inlet & outlet core flow distribution compatable with 
power distribution foreseen during the core lifetime. 

(ii) assure cooling of the control rod pins and hydraulic dashpot 
action (for slowing down control rod velocity at the end of 
its insertion) during all operating and accident conditions. 

(iii) DNB (Departure from Nucleate Boiling) should not occur 
under all operating conditions & transients ( i.e. Heat 
flux anywhere, to be kept below CHF (*3ritifal Heat Flux)«) 

(iv) Bulk boiling of the coolant should hot be allowed to occur 
in any channel. 

(v) The fuel centre-line temperature should not exceed 2500°C 
under any circumstance (fuel m.p. 2850 °C ) 

(vi) Hot-channel factors should be judiciously chosen; with 

the heat-flux increased by these factors, points (iii),(iv) 
and (v) above should be satisfied. 



i 
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4.2.2 CORE FLUX DISTRIBUTION 

The basic core flux distribution has been taken 
from the "Reactor 'Physics Design" section (Ch.5). Radial flux 
distribution has been assumed as a truncated Bessel's function; 
the axial distribution truncated cosine. Due to the high inbuiJ't' 
reactivity of the core,, control rods remain partially inserted in 
the core. At the BOC ( Beginning of Life/ of Core) the flux peaks 
as snown in Fig.(4.1) y (A). As the rods are gradually retracted 
due to burnup, the flux assumes a more sinussidal distribution. 
Finally at the EOL (End of Life)/ the distribution assumes profile 
(C). For calculational purposes, (B) is assumed as Sin(Z) 
and (A) and (C) as Z-Sin(Z) and ( ) . sin( .Tt-Z ) respectively. 

For comparitive purposes/ the normalized flux distributions have 
been shown as a computer printout. ( See p.27 of computer printouts) 

4.2.3 HOT CHANNEL FACTORS 

In the thermal design of a marine reactor core one 
is concerned with the most severe conditions in the core. If 
the design is adequate for these conditions/ it is certain to 
be satisfactory elsewhere. Hence the 'hot-channel ' is that coolant 
channel where heat-flux and enthalpy rise are at a maximum. 

(a) Nuclear hot channel factor ( Fq ) 

or Nuclear heat flux factor is defined as: 
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N Estimated maximum heat flux in core (based on nominal 
F U - fuel enrichment & dimensions), 

Mean Heat flux in core 

F^? = * F^ where, subscripts 'Z 1 & 'R' stand for axial 

& radial flux factors. This is the peaking of the flux in the 
axial Sc radial directions given as "leaking Factcrs'hin the CODE. 

F^ < 2.32 ; F^ < 1.57 for Sin(z) flux to F^< 1.82 

for Z-Sin( z ) and ( %-Z) gin.. (n-Z) fluxes. This causes 
1^^4.2224 ( We have obtained F^ = 3.4411 for Z-Sin(Z) flux.) 

This factor has already been included in the CODE. 

E 

(b) Engineering hot-channel factor (F ) . 

The maximum hdat flux may increase under the > 
following engineering uncertainties 

(i) Local fuel density being greater than nominal (factor 

E 

concerned - F^ ) . 

E 

(ii) Fuel dimensions being in excess of nominal (F d ) 

(iii) Local fuel enrichments being greater than nominal (F )• 

e 

f e s f e * f e * fI? 

Q *5 *d e 

For a PWR this gives a value of 1.08 to 1,10 
(Ref (26), pp. 170-177). But this approach is unnecessarily 
conservative and the effect of fabrication tolerances is evaluated 
statistically giving 

= 1.04 . ^Ref. (20)) 
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E 

(c) Engineering hot-channel factor for enthalpy rise )* 

This 'accounts for engineering uncertainties in 

the maximum enthalpy rise in the concerned hot-channel. 

E 

depends on the following factors 

(i) Fabrication tolerance subfactor, F^ accounts for fuel t 

enrichment variation, fuel weight variation & coolant flow , 
variation; it is evaluated statistically. 

(ii) Core Flow Sub-factor (E, ). accounts for the hot-channel 

An , Flow 

receiving less than its expected share of coolant flow. 

E 

(iii) Core mixing sub- factor (F, , . ) , accounts for flow from 

,Ah,*nixx; 

surrounding cooler channels mixing with flow in the hottest 
channel* - This is favourable countering the deleterious 
effect of the other subfactors. 

Typical values of above for a PWR are:- 

„E 


A-h Ah,F a b Ah, Flow Ah, mix 


1.08 


(d) Overall hot— channel Factor (F) 

E 

The engineering (F ) uncertainty factor would be 

F E = F E * F^ h ( 1.167 for PWR ') 

However, in the Italian ship Design (Ref(2) , p.598) F^= 1.22 has 

E 

been used instead of the above figure, giving F =l-.~269; this 
has been used in this analysis. This means if the heat flux 
is increased by this amount in the CODE, the core design steorrihd 



should sustain itself. (F^ has, as seated earlier, already bee 
included in the C' a ). 

The overall factor F, would hence be; 

N E 

F = F q * F = 4.357 

This implies the core should sustain a heat flux over its av-r 
value, of 4.367 times everywhere in the core. 

4 . 2.4 THE THERMAL HYDRAULIC C ODE 

The thermal design of a core is virtually 

impossible without the aid of a computer. In the computer C'—IE 

care has been exercised to entertain all thermal effects end 

variation of material properties with temperature. It lists 

output in a 3 -dimensional manner giving axial variations of 

temperatures & pressure for the rods at the core centre, ard 

then at discrete steps upto the core periphere. The analysis 

of the rod is for;- 

(i) Coolant integrity 
(iM.) Cladding integrity. 

(ii'O Coolant exit conditions and pressure drop. 

In addition for:- 

(iv) Radial flux variation in the core 

(v) Pumping requirements for primary system. 

(vi) Determination of core diameter (for 120 MWth power), 
and total number of rods. 

(vii) Temperature feedback for reactor neutronics calculations. 
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The importance of the CODE is felt when one has 
to do repeated analysis of a core with small variations in 
input, to arrive at an optimum solution. In addition/ cores 
of different reactors can easily be analysed since material 
properties are fixed. The CODE developed is versatile enough 
to be used either for stainless steel or Zirconium cladding, 
and either a square- arfayed • core or a hexagonal arrayed- 5 core. 

4.3 THERMAL AND HYDRAULIC DESIGN 

The fuel rod has been axially divided into 50 equal 
segments. The properties of materials in one segment are considered 
constant ana so also temperatures^, axially. The basic theory , 
governing the temperature distribution, heat transfer rates, 
pressure drops & energy balance is presented below, and done i 
for each of the 50 segments, successively. Temperature distribution 
in fuel Sc cladding are calculated using one-dimen$ ional radial s 
heat conduction. The pressure drops & heat transfer coefficients > 
are calculated using empirical relations for flow channels. 

In this core design, the square 

arrary has been used. For the 

fuel rods of diameter DCL an<^ 

pitch s , the flow ape a (A) pqr 

rod is given by j 

A = s ft *DCL^ (shadea portion^ 

. . . .Eqn. .{4.1 ) \ > 


4.3.1 FLOW CHANNEL 

i*~ DCLr* 



4 
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4.3.2 HEi.T GENERATION 

Fission heat generated is directly proportional to 
thermal flux (taking cross-sections as constant). The volumetric 
source strength ^’is 

111 T~‘ O 

“4 = 3.204 * 10" 11 * r 0 (r,z) . J/cm -s 

» *1 P 

Thus, q varies axially radially as: 

q''(z,r) =q_ * Cos ( Jt 2 ) * j ( 2.405 .r) Eqn (4.2) 

Hjt 0 R+d 

for za 0 to' H 

I I f 

r=0 to i-H^vhere q^ is the peak flux at centre of core, 

R= Core Radius 

d= distance from core face where 
flux =0. 

An important quantity required to be known beforehand is he at- flu: . 
generated per unit length of rod (kw/m ) q_^ . This is an average 

quantity for all rods. After a review of the reactors of the 
world (Ref. (14)), a figure of 11.0 KW/m was arrived at. 

This factor is related to the average volumetric heat flux thus- 

l 

q av = q av / % ^ where, r^ u = fuel pellet radius 

The volumetric power density at the core centre is •- ■. t 


f 1 I 



r N 

'R 


★ TTI^ * 

F z * 

,N 


i i t 

q 


where & F 

before. 


av 

,N 

Z 


are hot-channel factors as explained 



N3 JO 'Z 
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since by virtue of definition/ 


1 


it' 

o' 


1 f I 

all a 

o 

core 


Cos (JL Z ). 
He 


( 2.40 5 
R4*u ““ 


2 • rdr^qg 


where V= core volume 
we arrive at the following relations 


2 Jj££l + 2 (-S-) V^) 

i5, mR 


for a truncated Bessel ! s 
function flux-distribution^ 


JpJ 

— r *— rrr / Sin ( t ) , for truncated cosine 

/j li^ OT-T 

e 

1.8197 * % * H 


He (-ACos(A) -f B Cos, (B”) + Sin ( A) -Sin (B) ) 
/o~ truncated ( t-z).sin (ft-z) function 


, for 


where y A = H 
He 


and B= ( 


He-H 
2He ' 


Refer Appendix AS- I for the above inferences. 

4.3.3 TEMPERATURE DISTRIBUTIONS 
(l) In fuel rod - 

The basic heat conduction equation with heat generation 
is (in cylinder ical coordinates ) 

, 1 dT 


2 

d T 


dr 


dr 


1 I I 

.2__ 

K 

Fu 


= 0 
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with boundary conditions as 

T=To at r=o (fuel centre-lin^ temperature) 


dT 

dr 


= 0 at r=<o 


the solution of the above equation is 

I I I 


= To 


4K. 


Fu 


showing a parabolic temperature distribution with the maximum 
at the fuel rod centre. Hence, at the fuel outer radius (r F ^) , 

III 


T F = To 


q 


4 K. 


'Fu 


.... ( Eqn 4.3 ) 


Fu 


During the steady state condition, all of the energy generated 
must be transferred out to the cladding - since K changes with 
temperature, an iterative technique has been used to obtain the 
cladding thermal conductivity to within 1%. 


If K^T) = a+bT where a and b art constants, 

then we find that the value of Kp^ to be used is precisel 

To + T 

equal to K evaluated at the mean temperature F 


Fu 

i.e. in Eqn. 4.3, K pu (T) = K r ^( To+T F ) 


(Refer Appendix A4-2) 


(2) In gas gap 


The temperature drop across the gas gap is taken as 

( Eqn 4.4) 


Tj. -T = 
f g 


im r_ 

_2 ZJ*L 


2 h, 


gap 


1"1 2 
where gap = gap conductance = 0.6w r/r, m -K 

(p.356, ref (12)) 
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h is taken as constant throughout the core, 
gap 


( 3 ) In the claddinc 


Technique used here is similar to that in the fuel rod. 


The governing equation is ;• 


1_ (r dt ) = 0 ... Eqn . ( 4 . 4b ) 

r dr dr 


with initial conditions as 


T = T at r=sc^ (outer radius Of gap ) 

r) j>_3 


off = *■ r Fu 1 


r = r. 


where, K c ^ = thermal conductivity of cladding. 


Eqn. 4. 4b gives a simple logarithmic solution 

id 2 

T = T_. - i Eu , r_. 

g — z ~ - — — + In ( cl ) . 


( Eqn4 . 5 ) 


where r t outer cladding radius, 
cl 

Here again, K is a linear function of temperature and hence 
an iterative technique has been used to obtain to within 1% 


occur acy, 


K (T) = K ( T F + T S ) 

cl d — 2 


(4) In the coolant 


All heat energy generated by the fuel must be ultimately 
transferred to the coolant through the fluid layers near the 
cladding outer surface. Newton's law of cooling has been used 
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for this purpose 


/ 

heat transfer heat transfer ^heat transfer 
rate. coefficient area 



is a mass-weighted average temperature of the water in the 
flow channel. 

As, (r pu ) 2 A1 


A 


s 


2 % r el ^ 


thus. 



» > i 

r< r 


Fu' 


2- h -r 


c.l 


( Eqn 4.6 ) 


The method of determining (h j Is elucidated later in this 
section (4.3.5). 

Tq is to be known 


before we start temperature calculations. After that, 

T s ' T g' T f an <3 To are determined using equations 4 ,6, 4. 5 ,4. 4 and 4.3 

successively. Method of calculating Tg has been stated later 
(4.3.6). 

4.3.4 COOLAMT PRESSURE DROP 

We shall determine only the pressure drop in the core 
here. There are 2 types of losses occ ur i ng. 

(i) Frictional pressure drop along the channel 

(ii) entrance and exit permanent pressure losses, A Pg. 
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(i) A p f occurs due to the shear stress on the flowing fluid 
by the walls of tne flow channel. This is given by the relation 

^ P„ = f «■ * 2. V B ( Eq n 4.7 ) 

De 2g 

where, f-'=friction factor 

<il= segment axial length 

H = water density at Tg 

v B = average channel coolant velocity 
D e = Equivalent channel diamter, and 

/ 2 

g = gravitational accleration constant/ 9.81 m/s 

« ^ n = hrea = A_ 1 A 1 cUfined by 

ere ' e wetted Perimeter 2^r , ' 

c3 - (Eqn. 4.1) • 

The value of i depends on flow and surface conditions in 
channel. The flow conditions are characterized by Reynold's 
number / Re . 

Re = ^ V B D e ( Eqn ,4.8) 

P- 

where/ ji.= Coefficient of viscosity of water. 

Re should be greater than 2000 for turbulent flow that is 
desirable due; to the higher heat transport rates. In our case/ 

^ e = 237842.1. ( see p. 28 of computer printouts ). 

The second condition/ the surface condition is characterized ; 

by ratio of surface roughness (e ) height, to the equivalent diameter 
• Tne flow channel exhibits smooth tube behaviour if £ < t , 
the laminar sublayer thickness. Hence the cladding outer surface 
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smoothness should be so controlled that this criterion is met. 


The smooth circular tube friction factor ('F . ) 

cire 

is got by using tnis coorelation widely accepted in reactor 
work; 

£ cire = °* 184 * R e°' 20 (Eqn 4,9 } 

The Deissler & Taylor correlation for the correction 
of circular tube data to the rod-bundle geometry has been used:- 

f 2 , 

vr — = -3.475 + 8.053 (~) - 4.705 — ) + 0.9162<f_7 
r circ e De 

.... (Eqn. 4.10 ) 

This value of ’ has been used in Equation 4.7. 


(ii)<0>P E losses are the result of increased viscous energy 

dissipation resulting from increased turbulent activity 
at inlet and outlet. 


P E = K 


±i 


2g 


where K = resistance coefficent for sudden contraction 
from large to small flow diameter; K=0. 5 at core inlet, 

K=l.O at outlet. 


,4.3.5 HEAT TR/JSTSFER COEFFICIENT (h) 

A two-step procedure has been obtained to calculate 'h 1 

(i) 'h 1 for flow in circular tubes in found first. 

(ii) A correlation correction is applied then to obtain 'h ' 

for actual reactor coolant flow geometry. No ' direte-t- -&nalyt''ical 
methods for predicting 'h 1 for turbulent flow . . . 
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exist. However, by use of an analogy between heat & momentum 
transport, dependence of 'h 1 on Re cetn be predicted. Then, 
by applying the results of a heat transfer analysis for laminar 
flow over a flat plate, an estimate of the prandtl number, P f 
dependence on heat transfer is obtained. In this analysis, the 
Dittus-Boelter correlation has been used 

N Uc ir = 0.023 * R e 0,8 * Pr°* 4 (Eqn.4.11) , 


where N,, . = Nusselt number for circular tube 

u cxr c 

R e = Renault's number (given by Eqn. 4. 8) 

n C |A 
p r = 

K 


Now, Deissler & Taylor correlation for the actual coolant 
cnannel is done, like in Eqn.4.10 


Nu 

N U . 
crrc? 


= - 3.475 + 8.053 (~) - 4.705 (4j- ) + 0.9162 (~ ) 


% 


, .. h D 

since Nu = e 


(Eqn. 4.12) 


h = 


K 

Nu *K 
D 


(Eqn. 4.13 ) 


4.3.6 TEMPERATURE INCREASE BY THE FLOWING COOLANT 


Coolant temperature varies from a minimum at the 
inlet to a maximum at the core exit. The temperature in our 
analysis is assumed to very discretely from one segment to 
another. The bulk temperature T at any point (required in 

Jd 


\ 
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Section (4.3.3 ) is determined from an energy balance between 
the total energy added by heat transfer, and the energy rise 
of the coolant; i.e. 

ft. C pi T =3-^4, A1 , 

where, ih = mass flow rate of coolant 
^ 1= length of the segment 

At= increase .in Bulk temperature from one segment 
to another, ( T 32 ~ T bi^ 

Thus, T = T_, + n r fu 2 .Al (Eqn 4.14) 

m p 

for the first segment T 0 ^= inlet temperature and thus T^ is 
evaluated. For successive segments becomes T B1 and the 

calculation proceeds. 


4.3.7 SPECIAL CONSIDERATIONS 

Two important checks are required to ensure proper , 
coolant conditions:- i 


(l) No Bulk-Boiling occurs anywhere in the core; 


To ensure this T 0 T 


‘SAT 


, the coolant saturation 


temperature at the local pressure. The critical location is, f 
of course the core exit for the central rod where bulk temperature 
is maximum and pressure the minimum. 
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(2) No local boiling or sub-cooled boiling occurs near the 

surface of the fuel rod. Whenever T at a local point, 

this phenomenon may occur like during peak loads. This may 

result in Boiling burnout of the surface. To assess this, 

the local heat flux , is compared with the CHF, Critical Heat Flua 

(or DImB) at which state blanketing of the heat transfer surface 

by coalesced vapour bubbles is just imminent, ZJiy further flux 

increase requires heat transfer through this blanket. Due to 

the low conductivity of the vapour blanket, a catastrophic increase 

instr'face temperature would result, leading to melting of the 

cladding. 


Determination of the CHF has been done using the 


widely accepted Jens and Lottes correlation: 

* M 

CHF * C < T SAT -T B )0 ’ 22 ...(Eqn 4.15) 

where, C, M = pressure dependent const ant s JC see sec.4.3.8 ) 


Whenever sub-cooled flow boiling exists to any 
extent, Ne*/fcon's law of cooling for determining T c is invalid 
(as 'h 1 in it does not reflect the increased agitation in the 
boiling film due to the vapour bubbles). Jens & Lottes 
correlation for sub-cooled boiling gives (in f.p.s. system) 


1 ' CO 

= 1.9 ( q exp (- P/900) + T gAT 


where, P = local coolant pressure 

II I I I 

q ~ q - 


r , 2 
fu 


2 r 


cl 



S3 


The CHFR (CHF ratio) defined as ratio of CHF and q. B , should 
hence be grdater than 1 to satisfy condition (2) above. 

4.3.8 OTHER EMPIRICAL .CORRELATIONS 

The thermal conductivities, viscosity, density 
and specific heat at constant pressure are all functions ,of 
temperature. As temperature in the coolant channel continue 
varies, correlations for them would be required . 

They are (in c.g.s. system) 

K fu< T #) = 0.10308 - 0.14186 E-3 * T + 0.73263E-7 * T 2 , 

K (T) = 0.14194 + 0 . 6924E— 4 * T, 

* 

* \ 

for the coolant :- 

K (T) = 0.27656E— 2 + 0.29731E-4 * T - 0.10152 E-6 * T 2 . 

U (^)= 0.329E-2 - 0. 12479E— 4 *T + 0.15296E-7 * T 2 ; 

"S (T) = 0.7.38 + 0.6445 E-3 *T - 0.45827E-5 * T 2 ; 

C (T) = 17.677 - 0.1104 * T, + 0.2338 E-3 * T 2 ; 

P 

certain other phenomena depend on local pressures. They 
are correlated as under ( .P in p.s.i. ) 

T SAT ^ = 0*23357E3 + 0.6129E-1 * P - 0.5666E-5 # P 2 , 

C(P) = 0.11312E1 - 0.78E-3 * P + 0.3312E-6 *p 2 - 

0.56355 E— 10 * P 3 , 

M(P) = 0.48141E-1 + 0.22471 E-3 * P 4- 0.8302E-9 *P ? , 

(C 6c M as in eqn 4.15 ) 
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Finally/ a series expansion has been used to evaluate Bessel’s 

functions J (x) & J, (x) and modified Bessel’s functions I (x) S 
O 1 o 

I^(x) (from ref (40)). 

4.4 RESULTS AND COMMENTS 

Analysis of core performance Q f merchant ships 
SAVANNAH, MUTSU and LENIN was done using the Thermal Hydraulic 
CODE & results obtained were within 1°C accuracy. A pellet 
diameter and clad thickness of respectively 0.858 cm Sc Q.,66 mm 
was chosen after an optimization study using ref (14). Results 
for various combinations of pellet dia. and clad thickness for 
a constant volume ratio of fuel Sc water (for satisfying Reactor 
Physics considerations) are given below; the design chosen 
(case 1) can hence be justified as being optimum:- 


Parameter 

Case 1 

Case 2 

Case 3 

Case 4 

Pellet size 

0.858 

1.000 

1.000 

0.700 

Clad thickness 

0.066 

0.066 

0.100 

0.066 

Pitch 

1.30 

1.766 

1.810 

1.262 



— 

182.41 

186.99 

130.34 

Core Power Der 

r 

*63* 51 

34.42 

32.76 

67.41 

Max fuel Centr 
line Tempr. 



e~ 

1696.8 

1678.6 

j 

1688.4 

1725.3 

Min. CHFR 

obtained 

2.701 

3.334 

3.530 

2.313 

Max. burnout 

1.615 

2.946 

3.000 

1.809 

Burnup flux 

1469 8 

16 982 

16530 

10685 
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The thermal Hydraulic CODE can be used for various flux 
functions Sin( z ) , 2 .Sin (2 ) & ( it -z ).Sin(t -z ) by 

feeding an input ’flux code' 000,111 or 222 respectively. 
Some results obtained by this exercise ensue as follows 


Parameter 

Sin ( z ) 

z . Sin ( z ) 


Minimum CHFR 

3.248 

2.701 

2.797 

N 

Total F^ 

3.705 

4.367 

4.367 

Max. fuel C-Line 
tempr . ( 0 C ) 

1541.4 

1696.8 

1696.8 

Max. dadoing tempr. (°C) 

335.86 

335.83 

335.88 

Flux beyond wnich bulk 

boiling/burnout occurs 
(times of aver?ge flux) 

1.615 

1 ■ ■ 1 ' 

1.586 

1.586 

Coolant exit tempr. 
from hottest cnannel(°C) 

316.33 

316.92 

316.91 


Results obtained from the CODE are given in 


pages 28 to 36 of computer print-outs; flux distribution 
illustrated is Sin( z). For comporitive purposes, fuel-rod 
conditions for the hottest channel for Sin ( 2 ) Sc ( % -z) . Sin( tc 
distribution are tabulated in pages 46 to 48 of printouts. 


********* ******* 
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5.1 INTRODUCTION 


A nuclear reactor is an assembly of a proper 
size of fissionable material and other associated materials 
essential to sustain a chain reaction. The power generated is 
directly dependent on the neutron flux. It is of prominence to 
understand the behaviour of these neutrons in the reactor. 


In this chapter, we consider those reactor -core 
calculations that provide k^^, effective multiplication factor- 
and the stationary neutron-f lux distribution at Steady-state 
operation, k ^ is required to establish nuclear composition and 
configuration, which satisfies criticality and control requirements. 
The steady-state flux distribution must be known to calculate 
reaction rates and power distributions that are needed for the 
thermal & shielding design of the reactor. 

Evaluation of k ef £ necessitates the determination 
of the 4 factors 'n» e »p, f and non-leakage probabilities and 
* For sketching the flux distributions, 2-group 2-region theory 
has been used with water as the ref?Iecto:n Heterogeneous effects 
have been considered at all stages of calculation. 

5.1.1 Neutron Cycle : 

A typical history of 100 fast fission neutrons 


and how they sustain a chain reaction is illustrated in Fig (5.1) 

All important neutron interactions are indicated in it; the 
parameters used are applicable to the present design. 



! 

1 

t 


100 

FISSION 

NEUTRONS 



\V'=2.5#/fission 



\ 



FIG. 5-1 CORE NEUTRON CYCLE 


i 


I 
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A computer COPE has been generated for determining 

k and the flux distributions. The amount of 'poison' required 
err 

in control rods has also been estimated. All parameters used 
have been corrected for the operating temperature and pressure. 
The inputs like core dimensions and temperatures of fuel & 
moderator have been taken from Chapter 4 • An enrichment of, , 

m 

5 . 0 % has been chosen (chapter 3 ). The scheme of fuel loading , 

in reactor is as per chapter 3 . 


The scheme of computations is as follows. Input 
data are enrichment, fuel- rod & reactor dimensions, _ref lector,., 
thickness & average fuel & moderator temperatures. The only 
variable parameter fed in is the volume fraction (of moderator) 


of Pois on, P £ rac * Pf rac is changed by interpolation till the 
critical determinent = 0 ; this will denote a steady-state 


reactor operation with the nett control poison of p frac * The 
flux in both core and reflector are sketched for the above.,, 


5.2 CRITICALITY OF THE REACTOR 
5.2.1 Evaluation of 

Considering the reactor core as homogeneous will 
be err«ru eous as neutrons may undergo several collisions within 
a fuel rod before encountering the moderator. Calculation of 
kj J1 £=-^p e«pf still holds*; but differences with respect to . , . 
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us considerations should be noted. 


ETA ( 


This gives the ratio of neutrons produced to 


that absorbed in the fuel. rod. Since absorption cross-section 


of Oxygen is essentially zero, 

•V) ^25 ^£2 5 

T :Sr a25 + ^a28 


(Eqn 1) 


where, $25 = 2.47 and ^ 25 ' ^*£25 ^ fe P en< ^ ent on fuel enrichment 
(hereafter), and cross-sections corrected for average fuel 
temperature .obtained from chapter 4. 


THERMAL UTILIZATION (£] 


This gives the probability that a thermal neutron 
is absorbed in the fuel itself with respect ,to the who.le core . 
i.e. ratio of neutrons absorbed in fuel only, to that absorbed 
in fuel, moderator & other structural materials. 


f = 


Xfu ■ ^Tfu • v fu + 3 Ln- f’lta V m +3 | s s 


where, all fluxes and cross-sections are average values . Since 


0 = 0 
jTss ^Tm ' 


f = 


’ . V _ 

‘afu fu 


Jf . V _ + (2 - v - V )*_ 

*hfu fu am m Tiss ss *T 


. . . Eqn ( 2 : 


where = ^TjM* is called thermal Disadvantage factor (TDF) 


(ranges from 1.00 to 1.05). 
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This factor is peculiar of heterogeneous reactors/ taking 

into consideration reduction of thermal flux in *£u£.l- rod- 
due to self-absorption. As -ffis unknown, 

f is evolved through lattice functions for cylinderical fuel 
rods:- 


l_ =SaL_!sL_i£lgs * F(x) + E(y , 2 > 

f ^'af u* r u 

.... Eqn ( 3 ) 

where F Cx) and E(y,z) evolved from diffusion theory, are 
lattice functions defined by 

x * I (x) 

F (x) = — 


E (y,z) 


2 I^x) 
2 2 

_ - 

2y 


i (y) K i -(z> 


O 


+ *o (y> 


I^z) 


^(zTT^Ty)" - k 1 (z) i-^y) 


where , x= 


L ' 
fu 


y = 


m 


m 


a = fuel rod radius . 

b = radius of equivalent fuel cell = 
thermal diffusion lengths. - 


Pitch 

y a 


The ^ above equations, however, are valid for a<&- b and h ev,tfc - 
inaccurate for closely packed lattices. For x, y., 
individually being -C 0.75, the following series expansions 


qj-q quite accurate :- 

F(x) - ■■ ' 


2i 


= 1 * 4 ? 


r (?) 


I 2 


Eqn (4) 


E(y, Z ) = 1+4 i - 4—5 -In <5-> ’ 


2 (i-2) 


2 2 
2 -y 


i*2 


. . .Eqn. (5) 


, t 



67 


*£' is the most important factor that controls k ef £. Control 
rod positions greatly alter it. The above treatment does not 
include control poison. In Eqn. (3), the value of has to be 
modified to include the effect of Boron Poison. This is done 
by considering the poison as being equally dispersed in 
moderator and a term ( p frac * ^Lp) now 9 ets added/ where 
= macroscopic absorption (thermal) of poison & P £ rac is 
the volume fraction of moderator the poison occupies. P £ rac 
varies Tine^rly with control-rod position. 

RESONANCE .ESCAPE PROBABILITY ( p ) 

This is the probability that a fission 

28 

neutron will escape capture in resonances in U as it slows 


down to thermal energies. It is given by 

N -v -1-0 


p = exp .<— 

J m -%m m m 


) 


Eqn. (6) 


where 0- = fast fluxes 


Here we are neglecting fast absorptions in cladding material 

0 

and Poison which are negligibly small, m * above, called 

*eu 

the disadvantage factor for resonance neutrons is unknown. 

So p can't be evaluated through Eqn. (6). In a similar manner 
to f, p is evaluated through lattice functions. 

Defining '£' , resonance utilization factor 

f - Fast, absorptions in fuel resonance . ; . 

XT s 

Total fast absorptions \ 
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■ I ; x V_ * 
fu ■ eff fu fu 


l ‘e£f N fu V fu 0 fu + ""t ^kia V m ' 


(or) l/f r .= 1 + 


*V 

* Qm j 


$m .,_gn 


0 

kr-o 


N * I -V 0 
fu eff fu ^fu 


and from Eqn. (7) p = exp (- — ) 

1 - 1 


Eqn. (7) 

Eqn ( 8 ) 


As ;( 


0 . 


m 


w — - — ) is unknown, f can't be evaluated. We hence use 
*fu r 


the., following, lattice -functions (for cylinderical rods ) 
1 


f~ ^ N 
r 


fS * v 

l sm m 


fu ^eff V fu 


* F (x) + E(y;z) ..... Eqn. (9) 


where, F(x) and E(y,z) have same meaning as Eqns. (4) & (5) 

with x,y,z being functions of fast diffusion lengths here. ( 

Knowing f , p is calculated using Eqn ( 8 ). In eqn. (10) "f2 s 
xr sin » 

is taken directly from tabulated experimental values - for i 

water 1.46 cm ^ ( P.232 Ref. ( 23 ) ) and I eff the effective 

resonance integral, is well represented for cylinderical fuel; 

rods as 


Eqn. 

where A,C are constants (Ref. ibid) 
a = fuel rod radius 
•8 = fuel density. 


eff 


= A + 




( 10 ) 


( 10 ) 4 evaluates itself to 21.254 barns* 


i 


) 


\ 


i 
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t 28 . 


Note ; Though listed values are for pure U they are valid 

for slight fuel enrichments (upto 5%) as it does not 
significantly after the extent of resonance absorption 


in U 


28 


In Eqn. (9) the effect of Poison onlS should be 
• * sm 

considered*. The. effective values of "ILdff^ - shou 1 rl 


be taken. 


1 


(i. zr ) + 

L sm m 


(^s S p) 


ef f 


"S' + T> 

sm ~sp 




Eqn (11) 


and US j;, = ( ^ + 2 ) 

smeff sm sp 


FAST FISSION FACTOR ( e • ) : - 

All neutrons don’t succeed in escaping from fuel 
and subsequently slowing down in moderator. Some are absorbed 
in Resonances above fission threshold as well as upto the 
lower threshold of inelastic scattering in fuel# giving 
rise to fast, fissions# a favourable phenomenon. £ is the ratio 
of number of neutrons having first collision in moderator . 
to the total produced in thermal fission. For determination 
of § # the method of 3pinard# Fleishman & Soodak (SFS ) 
modified for closely-packed lattices has been 'used. The .... 
neutrons are classified into 3 energy-groups in fue.l with 
corresponding neutron escape probabilities without having a 
collision in fuel as P^#P ^ & P^ defined by Wigner’s approximate 
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i” + 1 
at 


for i*l to 3 Eqn (12) 


where 




of fuel rod 


Surface area . 

= 2 ■« r^ for cylinder leal- rods. 
2^. = Total Group i cross-section. 


The 2 fractions of fission, neutrons in the 2 energy groups, 
with the first above and the second below the fission threshold 
of 1 MeV, are given as = 0.561 and ~ 0.439. The 
group contains neutrons scattered from the first group into 
energies below fission threshold given by fractions:- 




3 d-j-Afoui g ii 


it 


22 


9 c XS 

Z 2 t 


. 8 = 
' 3 


Eqn (13) 


33 

ir 


3t 


Now, c is defined as the sum of 3 probability terms as :• 


e = 


If'l 


1 





1 T3 




h t ( l-U-P^ P x ) j 


y * 




i-(i-p 3 ) P 3 


Eqn (14) 
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The,;, cross-section used in Eqn. 13 and Eqn. 14 are taken directly 
from p.402-406, Ref (23). 


k. _ 
inf 


Having evaluated the above k^ n £ is given by 


k. _ s \*P * f 
inf T 


5.2.2 LEAKAGE AND *£££-*- 


S ince the material & geometric 


k ±n£ -l 

buckling have to be equal for a critical system B 2 


th 


taking fast leakage into consideration 
= k inf * eXp ( ~ B X ) 


eff • 1+ L 2 . B 2 

th 


Eqn ( 15 ) 


In steady-state reactor operation , k £f = 1.00 

5 * 2 -*- 3 CONTROL REQUIREMENT!- 

X 

The amount of Poison required is determined 
from the CODE by trial and error method for varying values 
of P frac such that the critical determinent ( refer Eqn. 63) 
is a minimum. This will be the steady-state value of Poison 
required to maintain criticality. The inbuilt reactivity 
of the core is; determined by putting Pf rac = 00 $nd setting the 
core at room temperatures. The boron Poison required is 
evaluated by choosing a slightly higher value so that the 
reactor remains subcritical even in the cold condition. 
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5. 2*4 TEMPERATURE AND PRESSURE EFFECTS 


(i) Thermal group constants 

(a) The averaged thermal absorption and fission 
cross-sections are strong functions of temperature. Taking 
the Maxwellian distribution, they vary as 


a (t) 



Tp yO . 5 

T 2200 


V T) 




* 


c _ 2200 
f 


Eqn . (16) 


Eqn (17) 


where, T^ = 295.1 k T = average fuel temperature 

g_^(T) , g f (T) are non-—* absorption and fission factors^ 

respectively for deviation from Maxwellian distribution. 

For materials other than fuel, g (T) = 1.0 V T 

a 

The cross-sections for all materials have been corrected as 
above using following average temperatures :- 

Moderator (in core ) ........ From Chapter-4 CODE. 

Water ( in reflector) 280°£ 

Fuel temperature From Chapter-4 CODE 

Other materials in core Average Coolant temper ature 

(b) The diffusion length for water is a strong function, of 
temper ature and pressure. An experimental value for 1*^ 
given by Deutsch ( p.92. Ref (22)) has been used* 
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*■ /« 2.719 

L.,(T y pr) = 

th *§ f(T) 


* ( 


T+273.1 


295.1 


0.5 

-5 


1.0614 


1+0.0614 ( — 2 9 5 7 i *' 1)Q * b) 


i0.5 

) 

Eqn (19) 


where ,. pr ■= water pressure, psi 
T = water temperature/ °C 

(<i>) = water density at p r and T 


(c) For water. 


2 °< N, number density, which in turn is 

O- 

<K ^(T) 5 o< vgfr) . 


As, iJ(T) = 0.93796 + 0.64448E-3*T - 0.458271E-3 *T 2 gm/cc 

Eqn. (20) 

where, T= water temperature (for 2000psi pressurised 


water) 


So, 2 (T) =*5(T) * (~) 0 * 5 *2 (To) 


“a 'W ' • * 

o 

Next, to evaluate diffusion constant ! D for water, we use 

* t t 

the_ familiar equation 


D = Lt, * y 
fh ^a 


2 Eqn (21) 

We see that the D calculated in the above manner is 


reasonably accurate. 


\ 
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(±±) Fast group constants 


(a) The neutron age, t , in water is not 
that much affected by temperature as pressure, ic ; 

As temperature & pressure increases, density decreases 
thus increasing neutron age in water. We arrive at : 


x (T,pr ) = 


0 


Eqn. (22 ) 


4 (T) 


where 


2 x 

= 27 cm * at atmospheric pressure 

= 1 gm/cc | and temperature of 22.1°C 


0 


(b) The fast neutron age in core is dependent on temperature 
in an indirect manner. A temperature increase results in a 
decrease of thus increasing the age. This has been elaborated 
in the subsection that follows. 


(iii) Heterogeneous Core constants: 

(a ) In a heterogeneous lattice, the 
average absorption cross-section 2 * is 

2 = 'V + ( y + 2 ) vl & .2 . - v 

as afu fu ap ; ™ ' ass ss 

.. Eqn (23) 


and 


2=2 - v + 2 +2 

*r trfu fu tr»»i 


V. + . ’ m ^“trss V ss 


- Eqn. (24) 


'trp being relatively very small may be neglected. 


* Hereafter subscript 1 s 1 ..stands for thermal/slow neutron 

parameters. Subscript *£' stands "for. Fast neutron parameters. 

1 v, 



(b) The extrapolation distance (E.D. ) of the core is then 


E.D. 


0.71 


•d. 


^tr 


Eqn (25) 


(c) Since average diffusion coefficient, D , c\f 'Loctice varies 

5 

inversely as J tr , it follows from Eqn. (24) that 



V 

D 


fu 

fu 


* 


V 


l*L 



+ 


V 

ss 


D 

ss 


(Dir) D 

o 


1 


fu 

, Yh 

v 

i ss 

fu 

d m 

*■ D 

SS 


Eqn ( 26 ) 


(a) Now Ls , the thermal diffusion length of core, is determined 


from Eqn (26) &_(23) 
D 


€ 


S as 


Eqn (27) 


(e) s '«The neutron age in the core is determined from experimental 
values for age of a Ujr^O mixture (p.206 Ref . (23) ). In the 
above case, a change in temperature ( & pressure ) has the 
effect of reducing the atom density of water relative to 
fuel which is equivalent to increasing the metal- to- water 
ratio in volume. 

As 'S is given by Eqn. 20, the ratio of U-f^O mixture, say, 

X = Vol.of U * 

Vol. of H 2 0 S(T) 

( 23 ) 

Using Fig 6.4 of ref ' , we arrive at the following correlation 
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X 1 = 19.8 + 12.5 * X - 5.0 * X 2 y for X X <^0.9 
= 27.0 + 6.27 * (X-0.9 ) j for X> 0.9 


(f) In Eqn (26) D is taken from Eqn (21) but D * D and D are 
^ m ss 8 9 fu 

still unknown. These are individually calculated as 


1 

D s = t (2j +S s ( iqr» 


Eqn (28) 


(g) Finally/ the fast diffusion constant/ in core since 

hydrosjen scattering cross-section varies significantly with neutron 

energy, an artifice is used for core transport cross-section. The:v». 

using N j 

V e££ __ Eqn (29) 

^af = ln( ^JL ) 

m 

and Eqn. (28), is calculated. All the above calculations are 
shown in Appendix A5-1. 


5.2.5 HETEROGENEOUS VERSUS HOMOGENEOUS :- 

We are now in a position to comprehend what 
effects would we. ihave incurred had a homogeneous treatment of 
the reactor been done. As ^>1 } f (hetero) ^ f (homo) due to 
the self- shielding effect of the fuel rods. However, p (hetero) 
p(homo ) as the above effect reduces the number of neutrons 
that are parasitically absorbed by fuel resonances. Thus, decrease 
in f is more than offset by increase in p . So, (f.p) hetero 
(fp) nomo. 
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The fission neutrons pass through a region of pure fuel where they 
are likely to induce fast fission, before emanating to the 
moderator. Hence, ( e) hetero ( s ) homo. This effect is not 
that pronounced as the other factors above. T does not, however, 
change. 

In view of the above inequalities, for 

heterogeneous lattices is greater than its homogeneous equivalfrt. 

A homogeneous treatment, though much simpler would consequently 
give erroneous results. 

3. FLUX DISTRIBUTION 

The thermal reactor consists of a heterogeneous 
multiplying core surrounded by a non-multiplying water reflector. 
Since the core consists of a closely-packed square lattice of 
fuel rods, it can be considered homogeneous for determining flux 
distributions in core Sc reflector. Despite its limitations, 2-group 
approx imatiom has been used, as it is the simplest that will give 
reasonably accurate results for multi-region calculations. 

The diffusion equations for fast Sc slow neutrons 

t%nju 

JtS" given by :- 

D f7 2 0 f - ^ rf ' 0 f + kln£ • 7 ac^s = 0 Eqn (30) 

h' 

D sF 2 - Xas-^s + P'^rfh = 


0 


.... Eqn (31) 
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where, all symbols have their usual meanings , 


2 as = thermal absorption cross-section 

= s lowing-aown (fast removal) cross-section. 

In core, 1 and in reflector k inf = 0.0 

The solution of Eqns (30), (31) are solutions of wave equations 

7 2 jzf f + B 2 = 0 
V 2 + B 2 0 s = 0 

This implies the coupling of the solutions by the condition the. 
coupling- coefficient , S , between fast & slow fluxes is 


S = 


0f Z,__ + D B 2 


0, 


~’as 


-v 


Eqn (32) 


rf 


The criticality condition should be satisfied; 
k 


i.e. k nj> = l = 
eif 


k. _ 
inf 


Eqn. (33) 


( i+sl ) d+-§-) 

x| 4 

where x = l/L 
s s 

x f “ 1/L f 

We shall use subscripts '1' for core and '2* for reflector. 
Eqn (33) defines 2 values for B^ refer f i9-5.2) thus 
,2 


2 • - - vT - 2 2 — - 2 


B 


- <x if + x i s > + 


x lf + x ls> + 4 x lf Z 1 


s 


Eqn. (3&) 
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B. 


( x 2 .,-* x 2 „) - {/( x 2 ^ + x 2 _) 2 + 4 -x n „ -x, 2 • (k^^-1) 


k if A ls‘ 


\1£- is- 


; le if • inf 
Eqn. (35; 


B- = - < 

2 2s 


’ 2 2 
B 2 “ X 2f 


.... Eqn . ( 36 ) 

Eqn. (37) 


Considering wave-equations are separable in both core & reflector/ 


0 


Cos (- 


% 


H + 2RS 


-) * F(r) Eqn. (38) 

D 


R_— R 

where/ RS = reflector saving* = (—=—- — )- ( 1-exp (— — 1 ) 


2s 


2 s 


F(r)- is the solution of the equation 


d F 


r + ii + < fi2 - < H^RS 


dr 


since B has 2 values in both core & reflector (given by Eqns(34) 
to (37) ) / to define F(r) completely, a total of 4 equations will 
have to be satisfied 

d 2 F lf + 1_ d F lf + 1 2 = 0 Eqn. (39) 


dr 


n 2_ 
d F 


If 


dr 


dr 


+ JL. — F lf - m 2 F _ = 0 Eqn (40) 

r dr 1 


12 
d F 


dr 


„ d F 

, .1 2f 

2 r dr 


2f 




F 2f = ° 


. Eqn (41) 


* From p.69, ref (24) 
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d F 2f -l- 1 aF 2£ - H 2 * F =0 Eqn (42) 

.2 r dr zt 

dr 


Identical 4 more equations are there for F 

1/ m, j* s and ^ are all +ve constants s~ 

2 2 
,„2 = , B ' 2 „ ( — 

1 V H+ 2RS ^ * 

2 2 

2 _ R * z , £ f 

. B 1 "" ( H +2PS ' 

9 9 71 ^ 

B 2 + * H + 2RS^ 

2 -t 2 71 2 

P 's B 2 + H +2RS 


and F„ as well. 
2s 

Eqn (43) 

Eqn (44) 

Eqn (45) 

Eqn. (46) 


The solutions of Eqns (39) to (42) are Bessel's functions. The 
complete solutions are as under 
In core 

F ls ( ^ = A” J (Ir) + C'l (mr) Eqn. (47) 

oo 

F f (r) = A-S^^dr) + CS^ I (mr) Eqn (48) 

In reflector 

F (r) = E*I Q ( ^r ) + G *k q ( V~ r ) + H-I Q (^ f r) +M.K Q (n f r) 

Eqn. (49) 

F 2f (r) = S 2 ’»( H * I Q ( P^r ) + M*k Q (P' f r)) 


Eqn. (50) 
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where s given by Eqn. (32 ) are 

,2 


s i ■ 


s i - 


2, + D • B\ 

las Is 1 

P ‘-^lrf 

^as +D 1S *i‘ 


xr~i 

P-JL 


Eqn (51) 
. Eqn (52) 


lrf 


s 

2 



= 0.0 


v n ^ 

^2a.s + 2s" 

Vi 

*-2rf 



Eqn (53) 
Eqn (54) 


Eqns. (47) to (50) contain 6 unknown constants A, C, E y G, H 
and M. These can be determined by the following boundary 
conditions 


F 2f ( V = 0 


F 2s (R 2> = 0 


’WV = F 2 f (R l } 


F ls ( V = F 2s (R 1> 


If 


d F 1f (R 1> 
dr 


= D 


2f 


d 

dr 




dr 


■ ° 2S 


dr 


Eqn (55) 

Eqn (56) 

Eqn (57) 

Eqn (58) 

.... Eqn (59) 


Eqn. (60) 
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Applying Eqns (55), (56) to (45), (50) gives 


E 


*0 ( ^ R 2 


0 


Eqn. (61) 


H 


M 


K q ( HjE R 2 1 

V 


W 


Eqn. (62) 


using Eqns (57) to (62) in (47) to (50), we arrive at the 


following 


A 

C 

*11 

Y12 

Y21 

Y22 

Y31 

Y32 

Y41 

Y42 


G M 


Y13 

Y14 


r— — 

0 

Y23 

Y24 


0 

Y33 

Y34 


1 0 , 

Y43 

Y44 


0 


— 


- -J 


.Eqn. (63) 


wnere, Y. . 's above are elements of the critical determinent 

*1 J 

as given in the CODE ( line nos 23900 to 25900 of main 
program \ , ui.. 2 1 di ^ 5 ~r r'.tl-c _l.j tii« onl* «.r t 

•.p- ■■ the critical deter mi r, --nt 

; only then will A ,C,G & M will have non- 
trivial values as the 4 equations are homogeneous. 

Next, the first 3 equations of Eqn (63) are taken 'and constants 
C, M and G evaluated with respect to A $sing Crammer's 'rule*, 


i . 


given by 


Eqn (63) =0 
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Finally, A is determined by the following method. Since the 
average thermal £hx in the core is given by the relation 


+H/2 


%R 2 H 


< A * Jq ( lr ) + C ' Jq (mr ) ^ • Co s ^ 


— H/2 0 


2 % rdrdz 


from A 3.1- 3 , this works out 



2( i> 


I 1 (mR 1 ) 

mR^ 



f _ JLiL ) ( 

^ (H+2RS )*2 t I 


Hence, 

A_ = 1 * S 

2 sav “i 118 !’ + 2( £) I l (mR l ) 2(H-2RS)Sin^ H 1 ‘ + f RS} ) 

A "«i 

Eqn (64) 


It should be recalled that ( ~ ) in Eqn. (64) is known and we 
obtain 'A 1 as the number of times of 0 sav « From this, the other 
5 constants can be evaluated and since all constants in 
Eqns (47 ) to (50) are now known the radial flux functions in 
core Sc reflector can be easily sketched and is shown in Fig. 5.3 




5.4 RESULTS 


A simple one group calculation of k^ nf versus number 


.28 


den. ratio of U & water as per p,306, ref (23)/ gives a maximum 
at a ratio of 0.156 for the nuclear parameters chosen/ this 
implies a pitch of 2.12 cm for a fuel rod diameter of 0.858 cm a) 

3 

a core power density of 24.52 MW/m . This is clearly undesirable 
for the compactness desired. However/ a pitch of 1.30cm chosen/ 
though resulting in lesser moderation is a compromise/ with the 
operating point falling only 8% from the peak. 

Core criticality has been achieved at ^^^=0.992945 
with the aid of a 10.99 cm thick water reflector. The various 
intermediate parameters/ criticality determinent & others are she 
in file CHECK-OUT ( p.39 of computer print-outs) . Other nuclear 
parameters like cross-sections/ etc. are displayed in file K.OUT 
(p. 37—38 of Print-outs). Fig 5.3 maps the radial flux variation; 
a rise in thermal flux at core edge due to moderation by the wate 
reflector is clearly seen. 

Results of reactor control calculations (refer chap 1 

I 

evolves a volumetric Poison fraction of 0.92148E-3 for a clean ; 
reactor with a maximum of 2.35E-3 (with all rods in) to drive 
the reactor subcritical by 5% Ak/k in the cold condition. File 1 
PSN.OUT (p.40 of print-outs) lists some parameters discussed | 
later in Chapter 7. ’ 


************* 
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REACTOR SHIELDING 

RESUME 
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6.1.1 Types of Shielding 

6.1.2 Shielding Requirement a 

6.2 SALIENT -FEATURES 

6.2.. 1 Shielding Design Characteristics 

6.2.2 Dose Rates Permissible 

6.2.3 Optimization Techniques 

6.3 RADIATION ATTENUATION IN SHIELD 

6.3. 1. Nautron attenuation calculations 

(i) Fast neutrons 

(ii) Thermal neutrons 

6.3.2 Gamma Ray Attenuation calculations 

(i) Core gammas 

(ii) Capture gammas 

6. 3.3 Emissions from reactor effluents 

6.3.4 Total Dose -rates 

6.3.5 Weights. -of. shielding materials 


6.4 


COMMENTS 



INTRODUCTION 


Shielding of a marine 


reaetor is efjfaxMgr eater 1 importance than "its landbaeed 
counterparts.*. Here# the peculiar space# weight and dose rate 
limitations require an optimization study concerning the 
shielding materials and the differential shielding thicknesses 
on the btasis of dose rate mapping. Reactor shielding is 
essential to i- 

(a) Not only protect the personnel working in close proximity 
to the reactor system but also to 

(b) Reduce radiation levels in the vicinity of instrumentation 
and other equipment that may be damaged by radiation; 

(c) Reduce the. induced radioactivity of reactor components 
that may need to be serviced during shutdown (like primary 
loop components) 

(d) Prevent excessive radiation induced heating of the 
pressure vessel (PV)« 

(e) Limit the neutron embrittlement cf the PV. 

' The marine reactor shielding design would require 

a careful study, of propagation of radiation through matter , 
(i.e. shield material)# system analysis of activity production 
and release in reactor effluents* It also requires other 
aspects of the ship. & reactor design such as general plant 
layout structural studies (not done here), Ihe aim is to 
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utilize the large amounts of materials needed for shielding 
also for structural purposes to optimize on weight. For 
analysis a ^MARINE REACTOR SHIELDING COE® has been developed 
which has been referred to throughout . 

6.1.1 TYPES OF SHIELDING 

There are 3 Key factors in shielding design- 
cost/ shield volume & shield weight. Any one of these may 
decide, the type of shielding required. There are three types 
of shielding 

(i) Unit Shield - that completely surrounds the reactor, like 

in land-based reactors. - 

(ii) Compartment Shield - that shields the reactor proper but 

not all its components like in 
mobile reactors. 

(iii) Shadow Shield- that shields only the crew and important 

instrumentation, like in aircraft or 
space reactors. 

'i 

6.1 .-2 REACTOR SHIELDING REQUIREMENTS 

The .control of radiation from a ship-board 
reactor presents, in particular, the 2 following problems 
of design 

-shielding from, neutron & gamma sources in a compact geometry, 

layc?^. 

-management of safe release of gaseous activity to satisfy 


\ 
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the requirements of rapid access to the reactor containment 
after shutdown. Among these# only the former has been 
studied to some depth. 

The scheme of reactor shielding revolves 
primarily around these 5 considerations 

(a) Slowing down the fast neutron® / } n‘). 

(b) Capture of slowed down on 1 . 

(c) Determining the source of the secondary gamma radiations 
formed as a result of various neutron interactions in 
shield materials. 

(d) Attenuation of all forms of gamma radiation. 

(e) Determining the total dose rate of radiation emerging 
from the shield materials. 

Let us examine the above^ points; for high 
energy neutrons# use of elements of moderate or high mass 
number is effective# as they are able to reduce their energies 
through inelastic scattering. Since elements of low mass 
numbers are the best moderators# in the .form. of wat^r is 
used as a shield constituent for slowing down intermediate ■ 
neutrons * { IMeY) „ Hence a combination of a moderately heavy. 

i 

element like steel with water is used. 

As a general rule# most of the f-rays 
originating from the core get absorbed in the shield region 
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nearer to the core Itself. It is the secondary ¥«.rays, 
however, ( emitted due to thermal neutron capture within 
the shield) along with a tiny fraction of the core gammas 
that control the shield thickness required. The neutron 
attenuating material to an extent also serve this purpose. 
However, as ][ -my relaxation lengths are inversely proportional 
to material density, additional material of higher, density 
(like Lead) is incorporated to attenuate the higher energy 
V- rays. 

A word about the 'thermal shield'. This is basically 
to safeguard the PV. High gamma flux from the core leads to 
very high heat generation rates in any structural elements 
placed there. If it were the PV. then the resulting large 
temperatures & thermal stresses when added to the . already 
present high mechanical (pressure) stresses, would bring the 
total stress on it (PV) to an unacceptably high level. Hence, 
a thermal shield consisting of layers of stainless steel 
and water is placed between the core & PV. The arrangement 
is as shown in Fig. 6.1. The steel has a high melting point 
Sc excellent thermal conductivity. Energy is absorbed by 
attenuating Y-radiation and inelastic scattering, of fast 
neutrons. This heat is immense & is removed by the, coolant 
flowing into the reactor core through these shield layers, 
thus contributing to the nett available energy. from core* 




£ 


Fig. 6,1 


Thermal shield in the PWP . 
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The core- V energy is attenuated from ^ 10 17 MeV/cm^-a 
9 2 

to 10 MeV/cm -a and the fast neutron flux from 

0.28 x 10 14 #/cm 2 -s to 0.16 x 10 9 # /cm 2 -s before 

impinging onto the PV. Thus high thermal stresses and neutron 
induced embrittlement ( by high-energy neutrons) is effectively 
reduced. 

6.2 SALIENT FEATURES 

6.2.1 SHIELDING DESIGN CHARACTERISTICS 

Ihe scheme used is as follows:- 

(a) Thermal shield: stainless steel (35 347) and water 

alternately used. 

(b) Pressure Vessel: made of carbon steel # as it is required 

to withstand high mechanical stresses. 

(c) Primary shielding: also- called Neutron Shield Tank(NST) 

surrounds the PV-Made. of carbon steel 
v - & water alternately# it has important 

mechanical support functions besides 
shielding# like supporting the entire 
PV#GRDM# Cooling ducts# part of 
primary piping# etc and also to resist 
shock load acting on the ship. The 
concentric layers of steel and water 
• makes an optimal shield for both 
neutrons & if -rays. Lead is not used 
here due to its lower melting point 
& poorer mechanical properties. 



(d) Secondary Shielding - It surrounds the containment vessel 

Sc the complete primary loop (steam 
generator, associated piping, primary 
, pump, etc.)* It attenuates to the allowed 

dose, radiations coming from coolant 
loops ( reactor effluents)- mainly 
»-rays from 0 X activation) and those 
escaping from primary shield. It's made 
up of one layer each of Lead Sc Paly ethane. 
Concrete due to its much larger bulk for 
the same effectiveness has not been used. 
Study reveals that use of Boron based 
alloys here do not substantially aid 
in reducing the shield weight/bulk. 

Table 6.2 shows the thicknesses of the 
above shield materials evolved from the 
shielding CODE. 

6.2.2 PERMISSIBLE DOSE RATES 

The occupational individual accummulated dose 

is limited by the well known equation D = &*(N-18) rems, where 

N is age of the person. This implies a maximum dose rate of 

2.5 mrem/hr. For personnel working in nuclear sites the* limitation 

however is l t 5 rem/year or 0.75 mrem/hr. This is the design 

consideration in our shielding calculations. However a total 
(29) 

dose of 12 rem under emergency conditions has been allowed 

with doses^excess of 5 rems being adjusted over the next 5 years. 
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TA3LE 6 ± 2 

OPTIMAL SHIELD DESIGNED 


SHIELDING TYPE 

MATERIAL 

THICKNESS 

(CMS) 

REFLECTOR 

WATER 

10.99 

THERMAL 

Stainless Steel 

(Core Baffle ) 

1.27 

SHIELD 

Water 

13.90 

- 

Stainless Steel 

3.00 


Water 

3.50 


Stainless Steel 

8.00 


Water 

10.25 

PRESSURE 

CARBON 


VESSEL 

STEEL 

12.00 

\ 


Water 

32.00 

PRIMARY 

Carbon Steel ; 

16.00 

SHIELDING 

Water 

10.00 


Carbon Steel 

19.00 


Water 

31.50 

EMPTY SPACE 

(for housing components 
of primary loop ) 

Air , 

200.00 

SECONDARY 

Lead 

15.10 

SHIELDING 

Po lyethene 

13.70 


CORE RADIUS • • 67.145 cms. 

TOTAL SHIELDING OUTER 30.67 ft. 

DIAMETER 
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6.2.3 OPTIMIZATION TECHNIQUES 

The search for the optimum shield consists in 

evaluating the weight of various possible shield designs giving 

(28) 

same dose rates at specified points. B. Chinaglia , has 
sketched weight variation W for the optimum configuration, as 
a function of neutron (D ) to total dose (D n + Ecf) i.e. 

D / (D^+ DY). The curve proves that the minimum weight is found 
for D /(D +DY ) = 0.10. This gives the lead & polyethene 
thicknesses as mentioned. 

The optimum primary shield with respect to weight 
is composed by laminations of heavy metal located at about 
the middle of NST. No strong difference between the use of iron 
only or iron and lead is found. Hence only carbon steel is used. 

6.3 RADIATION ATTENUATION IN SHIEND 

The radiatiun from core includes a , (3 
and V-rays, neutrons of various energies, fission fragments, 
protons & a -particles from (n,p) & (n,cc ) reactions and 
neutrinos, oraly neutrons & gamma rays are considered for 
shielding design as they are most penetrating. The primary 
and secondary radiations from the core have been briefly j summed 
up in Fig. (6-3). We shall consider neutrons &V— rays 
separately. 




FjQ.5.a PRIMARY & SECONDARY RADIATIONS ‘SCHEMATIC . 


^ o 


6.3.1 NEUTRON ATTENUATION CALCULATIONS 

The theory of neutron penetration through thick 
layers is different from ordinary theories of neutron diffusion 
used for criticality calculations where the interest is in the 
behaviour of average particles. But here, the interest is in 
those exceptional high-energy particles that do not undergo any 
interaction w throughout the suield materials. The average neutrons 
that do the normal things are not very important as they stand 
no chance of getting f^r into the shield. 

Neutrons are absorbed in the shield in 2 distinct 
stages. First, high energy neutrons are slowed down by elastic/ 
inelastic collisions; these slowed-down neutrons are then 
captured in shield materials easily as capture cross-secctions 
are high at low energies. We shall now consider these 2 stages 
separately. 

(i) FAST NEUTRONS 


Theory of attenuation 

These attenuate as per the following 3 forms 

(a) Elastic Scattering- Here, neutron K.E. is lost by collision 

with target nuclei. Dependent on scattering 

nuclei mass, it is the most important 

form of energy loss for light nucleg. like 

H 0 where a single collision 

may result in the removal of the neutron. 

From l.eV to lOKeV, this form is dominant. 
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(b) Shadow Scattering: This occurs at energies greater than 

5 to 10 MeV wherein the neutron acts as a 
wave and is diffracted about the target 
nucleus. It occur es only when neutron wave- 
length is lesser than nuclear radius. At 
large neutron energies/ its cross section 
approached f [ R , geometrical nuclear cross- 
section. 


(c) Inelastic Scattering: Its importance is that the incident 

neutron is strongly degraded since it loses 
the recoil energy of target nucleus plus 

, r* 

energy equivalence of inelastic V -rays 
emitted. This is very important for heavy 
nuclei where allowable energy levels are 
numerous and smaller/ciistances between ground 
state 6c first energy level. Its threshold 
varies from 6 MeV for 0^ to few tenths- 
of- air HeV for Uranium. At high neutron 
energies, its cross-section also approaches 
HR 2 . 


In marine reactor design, we are interested in 
minimum weight. Limiting values of both (b) and (c) above approach 
. Hence total cross-section (*% ) approaches 2^*R 2 at high , 
neutron energies. As nuclear radius R=1.5 x 10" 13 * A 1//3 where 
A is mass number of the nucleus, the total mass cross-section > 
becomes : 

T oo 2 

*— ‘T 6.022x10^ (2ftR 1 

= - 


0.0852 

7 , 1/3 


2 y 

cm /g 





which means light elements like give minimum shielding weights. 


justifying what was stated previously. 
Equations used 


Neutron flux follows the exponential attenuation law 

•^“*7 

(x) = . exp (-e,.x) 

where x = distance from source 
0q = source strength 

2 = effective removal macroscopic cross-section 
For 2 to 12 MeV neutrons 

a = 0,35 barns 

of material with mass no. A. This is used to evaluated of 
unlisted materials like steel (refer Appendix A6.1) 

Applying inverse square law for spherical attenuation 
and considering 0£(R) as flux at core edge, 

J^'(x) = 0.£ ( R ) . exp (-J.r ),R 2 /r^ for r)> R 
As more than one layer is present, the exponential term for 
'n' layers gets modified to 




where 


v 

i= cross section , 
t^= thickness, of layer i 


The attenuation prior to the P.V. wall occurs 
in the water and steel of thermal shield region. Experimental 
Data (p. 153,(26) ) for this region show that when build-up is 
considered, 
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point kernel can ue represented tv the sum of 2 exponential 
terms. For a large planar source of core, flux distribution 
in x cms of water, 0i^ (x) , is 

(x) = 0o ( A*E 1 (ax) + (l-A) i E 1 (cx ) 
where, 0o = inciaont flux 

E^( ) = incomplete error functions 

(refer App. A6.2 ) 

A = 0.89 

a = 0.129 , c =0.091 (constants.) 

The attenuation of fast neutrons thro ugh z cms. of steel 
followed by x cms of water is then well approximated by 


0 f (z,x) = 0 f (x) * E^ ( ^ S .Z ) 


Where 5= removal cross-sec for steel. Flux thus 


ss 


obtained is multiplied by 
spherical geometry. 


R 


to account for the 


(R+z+x) 


Method used 

The method used in SHIELDING CODE is 'Diffusion 
Removal Method 1 . Here, the basic assumptions are 

(1) The penetrating component of source neutrons consists of 
high energy neutrons that suffer only a small energy loss 
and uncollided neutrons. 

/ 

(2) Neutrons suffering large energy loss are regarued as being 
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removed from the beam. 

(3) The removed neutrons are degraded in energy according to 
diffusion theory and don't travel significantly from point 
of removal. 


The equations explained above have been used in 
the code to obtain the fast flux at every layer. The flux 
attenuation through the shield material is sketched in Fig. (6.4). 

(ii) THERMAL NEUTRONS 


The source of thermal neutrons are 


(a) thermal flux emanating from the core. 

(b) neutrons removed from the , ^ „ 

fast-flux beam. These removed 

neutrons, forming local sources can be adequately 
described by diffusion theory. This source intensity^ 
0 f (r) * H .In any layer a distance L from core 

centre , 


(L) * I?/ r 2 * exp (-S(r-L)} 

where ^(L) = fast flux at inner face of the layer 
t^nd r y L. 

These neutrons can then be introduced into the highest energy 

group of an appropriate set of multigroup diffusion equations 

\ 

to calculate the lower energy neutron distributions. Here, < 


only 2 group . . . 



103 


theory has been assumed. Hence# all above neutrons are assumed 
thermalized. (This only results in slightly higher dose rates). 

If 0 (r) is the slow flux at any radius' r, 
s 


D-y^ 0 (r) - T. * 0 (r) ■+ 0^ (r) *2 = 0 

rX C as Q +* 


1 ■ l \ 

leakage term absorption .local source term, 
term 


where D = diffusion constant of layer 
s 

X - thermal absorption cross-section, 
as ^ 


(or) V 2 ,0 (r) - c 2 0 (r) + Q e ~^ r - Ij ) 

v S S 1 


o 


where 


, d 2 = N‘ 


Z /D 

i a.s' s 


Q = 0 F (Lh.L? 


D 


s 


Evaluating for y 2 = ~2*~dr ^ r2 — ^ (spherical coordinates# 

^ considering only radial 
variation) 


0 s (r) 


Ae 


-cr 


Qe 


-El 


2 cr 2 cr 


<c-S). Bjty^or 


( refer Appendix A6~3) 
where, x = (c-l£)r , y= x/ (l+2c/{c-J:)) 


As 


-cr 


cr 


<< e we modify as 
-cr 


0 s (r) = — ' 


Ae 


2cr 


+ O e ^ (c-X ) . E (y) e Gr 
2cr 1 


. . . eqn ( 1 ) 



This equation is difficult to evaluate. As the effect of Inverse 
square law is much lesser than exponential attenuation within 
a layer , the distance attenuation is neglected within the layer 
for 0 (r) evaluation. This is a resonable approximation If layer 
thickness <X-L. 


Now diffusion equation redqses to 
y 2 0 (r) - c 2 0 (r) + Qe"^ r = 0 

S 3 

where, Q = H e . ^ ^ L ^ 

D s £ 

The solution of this equation (refer Appendix A6.3) 
is 


*V r) 


Ae“ cr • - r 

+ ’ # Qe 

r “r 


( 


c 2 -£ 2 


(c 2 -X 2 ) 2 


which is easily solvable; flux thus obtained is multiplied by 
2 2 

L /r for distance attenuation. It must be mentioned that this 
method gives a slightly higher value for 0 (r ) than Equation (l). 

The thermal flux within every layer is thus 
obtained & has been mapped in Fig. 6.4 . These flux values will 
be required later to evaluate the capture gamma-ray sources 
(local) within the snield material. 
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6.3.2 GAMMA RAY ATTENUATION CALCULATIONS 
THEORY OF ATTENUATION s 

Depending on the )/- energy and absorber mass- 
number, the 3 absorption processes of interest are photoelectric 
effect, Compton scattering and pair production ( in MeV region). 
Secondary effects such as induced reactions, coherent scattering 
and Raman scattering have not been considered as their effects 
are negligible. For illustration, V -ray cross-section is 
shown as a function of energy for Iron in Fig. 6.5. 

(a) Photoelectric effect : Most important absorption at low 

energies «1 MeV) Energy is transferred to orbital electrons 

through Einsteins equation: 

2 

E = hf = W + ig mv 

\ \ \ 

Total energy work K-S- of elctron 
function 

It is a true absorption process as the V-ray gives up all its 

energy in one interaction. The cross-section for this interaction 

is proportional to z n / E 3 where Z= Atomic number, E-^ray energy 

in MeV. 

& n =• 3 to 5 depending on ray energy. 

(b) Compton Scattering : Occurs at high energies. Here, y -rays 
undergo billiard-ball type collision with electrons, 
according to compton's equation: 

E= Eo 

1+ Eo ( 1-Cos ©) 

0.511 


where © = scattering angle 
Eq= incident)f- energy 
E = Scattered\fenergy 




GAMMA ENERGY (MeV) 


FIG. 6-5* GAMMA - RAY CROSS “SECTION FOR 
IRON.* 


Fig .5.3 y ref (26) 
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It is nor a pure absorption process as fractional energy lost 
depends on incident angle. Its cross-section is proportional 
to Z/E ( nearly constant for all elements (except H 2 ) where 
Z/E =0.5;) . 

(c) Pair production : At higher V - energies, this star tiling 
phenomenon occurs wherein under the effect of the coulomb field 
of the nuclues, a high - energy (> 1.022 MeV) can give rise 
to an electron pair-positr on & electron, with any excess energy 
appearing as particle kinetic energy. Positron always gets 
annnilated by combination with an electron giving 2,0.511 MeV 
V-rays. Its cross-section is proportional to Z*(E-1.022). 

It is a true absorption process as annhilation radiation is soft, 
emitted isotropically and attenuated easily. 

The above 3 process are all added together to 
give total IT -macroscopic cross-section . For materials like 
stainless steel containing more than one nuclear species, the 
individual ^' s are added together to give the total. 

This is snown in Appendix 6.1(b). From the table of mass- 
absorption cross-sections for tf's (p-762, (21)), it is clear that 
except H , the most efficient y- absorbers are heavy elements 
like lead ( & iron, to a lesser extent) which have high cross- 
sections for both photoelectric absorption and p&ir production. 
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(i/) CORE GAMM AS 

The reactor core exhibits a complex V" ray spectrum. 
For shielding calculations it can be summed up as follows s- 


TABLE 6-3 


V -RAY SOURCE 

) PHOTON 

t 

ENERGY 

- (MeV) 


X 

l 1 

2 

4 


8 

Prompt fission MeV/Fission 

i 

t 3.45 
I 

3.09 

1.04 

0.26 

- 

Fission product MeV/Fission 

S.16 

^ 

1.74 

L— 

0.32 



- 

Total 

I S - 61 

cn 

00 

• 

1.36 

0.26 

! 

- 


Besides these/ capture V" -rays are also being continually emitted 
through thermal neutron capture in water_ cladding material/etc. 
All this is summed up as in Apprendix A6.4. The core has been 
considered as a large uniformly distributed source of volumetric 
density (V D ) dictated by thermal hydraulic design. As 1 watt 
power is S.lxlO 1 ^ f issions/cm^-s considering 200 MeV energy 
emitted per fission, the volumetric source strength would be 

i 

i in r. 

S v = 3.1 x 10 - V D * Bf MeV/cm -s 

where, Y\ 1 s are chosen energy intervals from 
1 to 8 MeV . 

since the core absorbs y-rays emitted by it, this self- absorption 
phenomenon has betn considered. This reduces resultant radiation 
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level in the core outer surface by F . times, where F . is 

pi pi 

self-absorption function at the if -energy E^ . 

„ Radiation intensity with self- ab sorp t ion 

£p _ 

Radiation intensity without self- absorption 
For a cylindrical source viewed from the side. 


F P± = 


% 


R Xi 


where, R = core radius 

=V -absorption cross-section at 
\ -energy E»A L 


If dose rates are to be considered at distances^ R, then the 
reactor can be considered as a point isotropic source of strength 

Vi 

= Sy at V fk Fpj. MeV/j t y where V = core volume, for 

V -energy E^-j_. This is the source strength used in the CODE. 

In the CODE, the point kernel method has been 
employed where the reactor is considered as a source as above. 
Since V-rays follow the exponential attenuation law, in 
spherical coordinates, the /-flux will be 


0 (r) = S i • exp (-2*) 


where r= distance from source centre. 


4 t 

For attenuation through different layers of thickness t^ and 
V* absorption cross sections 2^, the exponential term becomes 


- b 


“ 2 ( Vi } * 

all i 
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From the main beam of ‘ s , the photons being scattered 
out - of any particular narrow beam are replaced to some extent 
by particles being scattered- in from .adjucent beams. This 
results in an increase in partiale flux- density depending on 
V-energy & material and "b". Hence a factor , isotropic 

point-source dose-rate build-up factor (>1) is included in the 
f 1 uxequat ion . 

0 (r,Ey. ) = S. fpE ( exp(-b(a/ i ) ) 

• ^ 2 

4® r 

when the total dose is required to be calculated , the effect 
of every By. is to be summed up. This equation has been used to 
obtain the effect of various energy fluxes at different layers. 

(ii) capture gammas 

These rays emitted throughout the shield, are directly 
dependent on the thermal neutron fluxes. From this, and a 
knowledge of material densities & cross-section, the secondary 
gamma source-rate may be calculated thus:- 

Source (i, r) = * 0s (r) photons/cm^-s 

aS 

at layer i, a distance r from core centre. In most cases, a 
spectrum of ][-. rays are emitted. These have to separately be 


considered thus. 
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Source (1, r, S'-±) = ^ as (i/E^/i) * £f s (r) 

Here, again, E^is assumed to vary in discrete steps 
from 1 to 10 MeV. 

Since V -rays are generated throughout the shield, 
consideration must be given to the varying shield thickness 
available to them. So, each layer of the snield has been divided 
into a series of concentric rings of thickness > 1.0 cm and each 
of these are Considered as separate sources. Tne total radiation, 
outside the shield is obtained by summing the radiation obtained 
from each of these rings from all layers, plus that obtained from 
the core itself, following the Point-Kernel technique (described 
uncier "Core Gammas") of ][—. ray attenuation, for each energy 
interval, Ey^. This procedure has been adopted in the CODE. 

6.3.3 EMISSIONS FROM REACTOR EFFLUENTS 

The use of light water as coolant poses far 
lesser problem than other liquid coolant in this field. The major 
reaction of interest is t 

O"*" 6 (n, p) N 17 activation cross section 2.0E-5 barns * 

17 

N emits gamma rays of average energy 6.2MeV-The total activation 
energy (a) evaluates ass- 

= 8.2925 E-6 * 0^ Vj^^O/eff * exp (—0.09367 ♦ 'Z ) MeV/s 
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where , 0f v = virgin fission neutron flux 



effective water volume to be considered for 
shielding purposes. 


x = effective time constant for transport. 


In the abscence of Vj-^O/cff and J and shape & df mansions of all 
primary components (extra-core) required to calculate equivalent 
effective flux for snielding, this emission has ret been 
considered for reactor shielding here. However, since its magnitude 
is small, it does not introduce great errors in dose 

calculations. 


6.3.4 TOTAL DOSE RATES 

Once the neutron and gamma fluxes are known at 
all interfaces, it is easy to evaluate the total effective dose 
rates. For this, the flux for unit tissue dose (l mrem/hr) is 


used as simple conversion factors (CF) . Taken from P.7-66 of 
Ref (21) they are :- 

TaMe. F>-J^ Fluxes fo v unit tissue. 4-ose vates * 

- — — 

Particle £ flux for 1 mrem/hr dose. 

j 

I 2 

(a) Thermal neutrons $ 1500 #/cm -s 

, < 2 

(b) Fast neutrons J 7. 333# /cm -s 

(1 to 10 MeV) j 

i 

(c) V-rays jj 

577 MeV/cm 2 -s 



1 MeV 



RADIUS (cm) 

FIG. 6-4(1) NEUTRON FLUX DISTRIBUTION 
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Particle I Flux for 1 mrem/hr close | 

L L 


2 

MeV 

> 68! 

MeV/cm^-s 

4 

MeV 

jj 841 

2 

MeV/cm -s 

6 

MeV 

1 946 

2 

MeV/cm -s 

8 

MeV 

1 

j 1022 

2 

MeV/cm -s 


using these CF^, the total aose rates at different layers 
have been calculated & listed in the CODE results. 

6.3.5 WEIGHTS OF SHIELDING MATERIALS 

In a marine reactor the shielaing should weigh 
minimum possible while still being effectivejonce the shield 
dimensions have been arrived at after having analysed many 
combinations of varying thicknesses in the CODE, the approximate 
shielding weight has been calculated using simple geometry ; 

techniques. For ..the varying heights of different shield -layers, 
varying height-scaling factors have been arrived at, after 
extensive study of a comparative nature done on various marine ; 

I 

reactors (information of only tho-se 1 available). i 

The overall shielding weight (including water- | 

snields ) comes to approximately 838.2 tons. The weights of the f 

I 

primary system, Core, Control Rod drive Mechanism (CRDM) and coolant &; 

i 

piping are to be separately considered for the total reactor | 
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weight. An uncertainty error of 10% in the above determined 
value is anticipated' ■ . 

6.4 COMMENTS 

It has been demonstrated from the CODE output that 

the maximum dose rate from the secondary snield is o.735 mrem/hr. 

3 

and that from the primary shield 4.331 x 10mrem/hr. Overall 
shield outer diameter is 30.62 Ft and weight 838.2 tons/ a 
200 cms space for the primary loop machinery has been assumed. 

With an increase in this space, the shield weight will increase 
substantially as the quantity of lead required increases. 

The neutron flux at the core edge has been chosen in 
a conservative manner to ensure proper shield design under all 
operating conditions. Fluxes obtained from chapter-5 have been i 
incremented by 2.304 times (1.269 Engg. hot-channel factor; , 
1.178 control-rod peaking factor- Z Sin(Z) over Sin(z); 1.5408 , 

flux increase due to core burnup). I 

The use of point-kernel technique with build-up 

factors, however, does not fully account the fact that in-scattere<^ 

V"-rays is of a lower energy than incident ][ -rays, E^.. The scattered 

component consists of a spectrum of energies from 0 to Ey. Tni^ 

changes with distance traversed through shield. Further, the actual 

scattering process is anisotropic. So the dose rates obtained by 
slightly 

our analysis are/higher than the actual rates. As V-ray attenuation 

is correctly described by the Boltzmann's transport equations, 
improved analysis result-s when transport codes are used. 


***************** 
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In the aosign of a mar in- reactor/ the control 
ana instrumentation system is a key factor in provicing an 
efficient safe ana reliable plant. ~ warsnip nas stringent 
manoeuvring requirements which imposes large load variations in 
relatively snorter times. Consequently, trie system snoula be 
stable unaer all conaiticns of operation with no unallowable 
over snoots, n jvnowlouge of core reactivity balance ana temper atu 
coefticients is requirea to assess tne above. 

Tne reactor Physics CODE and subroutine Abi\GN 
in specific (p. 16 of computer jorint-outs) nave been referred to 
tnrougiiout this cnapter. 

7*1 K>« ACTIVITY CQfiPUTi hjenS 

Two important fission products Xenon & Samarium 
absorb tnermal neutrons substantially. In aauition the long-llvec 
fission proaucts also absorb tnermal neutrons. Thus, if the read 
were to remain critical it snoula nave sufficient reactivity to 
overriae tne a.cove aelcterious effects. In general, the only ; 
significant efiects tnese poisons nave on k &£f is a reauctlon of 

1 . . . I 

tnermal utilization, f to a lower value f . Tne poison reactivot 

, is 
P 

% = V'"' = -^ap * f /S afu 


. . . . Eqn. (7.1) 
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The variation of above parameters with temperature 


should be considered. 


135 149 _ © 

7.1.1 Xe and Sin Poisoning : 


Tnis is the most important fission product having 
a neutron adsorption cross-s^ctiun of 2.7 Mb(2200 m/s). T he decay 
cnain for tnis isotope is 


135 

Te 

/s 


< 30; 




6, 7hr 


. Xe cs 


Direct Fis&ion Yielu=0.061 Direct Fission 


p 2.6j§6 yr 


( Y I ) 


Yiel&*Q1003 

<V • 


(stable) 


Detailed explanations of differential equations involving the 
above cnain are given in many boo.cs &. will not be attempted here, 
The following results, nowevd , ensure:- 


- q i q- 

q + 


£qn (7.2) 


wnere, / x = X , a temperature dependent function 
X,. = Average fuel fit. sion cross-section 


= Average core thermal flux 


= Steady-state Xenon Concentration 


From Ref (23) 
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135 

Eqn. (7.2) substituted in Eqn. (7.1) (forZl^) gives the Xe 
activity at steady state. After reactor shutdown, Xe activity 
rises above this value. The maximum occurs at a time, t 


max 


max 




m 


( 


* A. 


X 


i + ( y i + V 

Eqn. (7.3) 


, = X 

I I 


where, = y \ /a 


ax 


The Xe concentration is tnen given by substituting t-t maX in the 
following Post-shutdown equation 


” t 

x ( t ) = X Q e X 




x i “ X x 


( o 


-Xt 


< Y I + V * Z f 


where, X Q - ^ - - 


ax * 


e 1 ) . . Eqn. (7.4) 


, Steady- state Xe 
concn. 


Y x r f • 


*0 = 


X 


, Steady- state I 


135 


concn 


obtaining 2 from Eqn. (7.4) and substituting it in Eqn. (7.1) 

3X 

the maximum Xe activity is evaluated. 

149 




Though a lesser poison, this has a thermal absorption 


cross-section of 40.8 Kb. It 's decay chain is 

Nd 149 P ^ Pm 149 P 


2 .hr 




54 hr 


149 

Sm ( stable ) 


direct fission yield = 0.0113 
(Y p > 
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The Steady-state Sm absorption cross-section, * 

is given by :- 


=Y p *K aqn - (7 ' 5) 

1^9 

Unli.ce Xe tuis is inaepenaont of neutron flux. The max. Sm 
concentration after snutaown is given by 


S 

max 


S + 
<y 


P 

o 


Xqn (7.6) 


occur in y as t- — ^inf, wnere 

S o~ Y p ^*f ^^aS l subscrijjt 'p 1 stands 

- „ I 

P q = Y^Sf 0 T /Xp ) for Prometnium. 

Obtaining £ g from &qn. (7.6) anu substituting it in bqn(7.1), 
the maximum Sm activity is evaluated. 

7.1.2 LGhG-LIVtD FISSION i-uOOUCTS 

Tnese also absorb tuermal neutrons. Many products 

naving varying cross sections for neutron absorption ana decay- 

25 

-times ensue by U burnup- (Refer Table 9. 2, Kef (24) )j; for 
computing tne long-term effects of trie accumulated products, 
an average cross-section of <±5 b per fission is chosen 

(p. 106, kef (24) ) . Tnis has been obtaineo neglecting the rapidly 

V. 

satin: ating products.- Since the. concentrate, n of 'these products ‘ 
i§ essentially tne burnt atoms of it’s macrpscopic 'cross-' 
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section/ S a £ , is given by 

= 3 frac * * 43£ ~ 24 E « n - <7 - 7 > 

25 

wnerfe/ B fr-c = Fract: *- on °£ u burnt at End of core life (see 7.1.4) 

n 2 ^ = number Density of U 25 ,^/cm 2 -s 

Substituting Eqn (7.7 in Eqn. (7.1)/ the activity of the above 
prooucts is obtained. 

7.1.3 TxmiF.no .TUko CCEr xiCIEi-TS OF K-oiTI/ITY 

Tne temperature coefficients of reactivity ( are the basic 

parameters required, for reactor St anility Analysis. As * e ff 

changes with temperatures of operation, a knowledge of oQ ' 

i. s 

at various reactor operating ivels both for tne fuel & moderator 
is essential. Unlive land-basea reactors, tne load variations 
here are from full-load to below 20% for automatic control. 
Consequently, ^ T ’s at various loads would be difxerent and require 
separate calculations. For convenience only full load, 3/4, h 
and ^ loads, as well as full-load at EOL (End of Life) have 
been considered. 

The basic equation for°^ T is 

'V'W - <* ln£ > + V p f ) +oC t (P t > ^ n - < 7 - 8 

c ^ T (i) = 1_. * di. 

i oT 


where , 


* (7*9) 
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i = Parameter involved 

Pp } = Fast, Thermal Non- lea,; age probabilities. 


Eqn ( 7 . cs ) can be expanded as below 

°S?^ K inf^ = V + T^ + T (p) + 

C< ' T (P T ) = . „T (cC - (l ) +oC ( 3 2 ) ) 

1+3 L t 1 

°^ T (P p ) = - 3 2 T t (pC t ( T t ) +<^ t (3 2 ) ) 


Eqn. (7.10) 
Eqn (7.11) 
Eqn (7.12) 


Eacn of tlie terms in Eqn. (7.10) to (7.12) can furtner be 
€ixpanaed but not aone here. (Rarer Cnap. 13, Ref (23)). 

Tne basic assumption for computations here is 
temperatures calculatea are uniform for the material? involved, 
throughout the core. In tne Reactor Physics CODE, every parameter 
is corrected for the operating temperature (except the core 
dimension variation waica is negligibly small), nence is 
determined by giving a small perturbation to tne operatiny 
temperature ana determining tne consequent variation in k ^ 
and nence°Sp using Eqn (7.9). Tnis was aone for temperature 
perturbations of botn fuel and moderator operating -temperatures 
inuiviaually , giving 4 ^ and oC. , the respective temper cture 
coefficients. The average fuel ana moderator temperatures of 
different loads were obtained from the Thermal Hydraulics CODE 
(chap 4) by varying the linear power rating of fuel rod 
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appropriately. Results obtained are summed up in Taole 7.1. 


Taole 7.1 Reactivity .Coefficients of Fuel & Moderator. 


1 LOAD I 

^frac 

| 

A - r_ rhV , FU£ 

eff , TEmPR 

1 

Full ! 

1 

; 

0. 9197E-3 ! 

! 

0 . 9929<* 

508.07 

3/4 | 

1 

1 

1 . 053F— 3 

0.39235 

442.27 

1/2 | 
1 

l.lti65iL— 3 | 

j 

0.33296 

3&2* 97 

J 

1/4 i 

1 

, 1 * 3 214 £—3 J 

0.99308 

I 1 

329.43 

; 

! 

ibUXj, | 
FUJLjL 1 

* it*— *-*-4. jr- 

1 
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.7.JL . 4 3uri i up_ & Reactivity B al an ce 

Knowleuye of tne fuel fraction burnt is necessary 

to obtain tile core burnup and reactivity balance. Tnis nas been 

acaievea using a term ('BURNT 1 in CODE ) , ie that fraction 

25 

of burnt U ? below wnich the core is tu.&ble to override tne 

Xe-Sm ana fissian proauct activities and thus becomes subcritical 

with all control rous removed (P_ =: 0.0 ). Tnis state of core 

frac 

is termed EOL (Ena of Life). 

Tne core inouilt reactivity/ reactivities required 
to come up from cold to not-zero-power^, ana zero to full- power 
are all determiuea by choosing tne appropiate temperatures of 



1^0 


fuel moatrator in tne CoDH and evaluating ,c . «.t the ZjCL , 

t^II 

core burnup is given by the relation:- 


Burnup - (wt. of U p*r ton ot U0 2 * 3 &ac ) 

(mwa/ton) , , , , T .25 , . _ . „ ... . 

Wt. of U burnt for producing 1 Mwa 


. . hqn ( 7 . 13 ) 


ana, finally, tee core endurance (in days of full power ) is:- 


bnaurance 


3urnu£i_* r JTotal wt. __p_f__UO 
operating power level 


. . . . iqn (7.1a) 


(Ketor line nos. 6900 to 7400 of p.17 of computer printouts for 
jc,qns(7.13 ) ana (7.14) ). 


Having known tne above results, tne reactivity 

balance of tne core is rnaae as in Table 7.2. Tne amount of HoISGiSi 

required in tile control rods snoula be suen that tne core remains 

subdr itical even in the Cola conaitien. H^nce, a P =2.35 h-3 

frac 

was found adequate to drive the core subcritical by over 4 % 

Ak (:c ^^=0. 9532) . Tnis means as tnere are 6 control rods per 
assembly , by a simple geometric computation, the rod diameter^ 
is 0.51 cm ( J .494 cm B^C pellet diameter), assuming a stainless 
steel claaaing of tuiacness 0.53 mm * 
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TV-wLx, 7. 2 CORE RE^CTI;>/_ITY _ _ Ahkti^ ! 

—A 

Cola to '0 1 jjowtr 9.6360 

■ O' to full power 5.4634 

j A.e , Sm poisoning 2.8563 

i 

•A-fc/Sm overr iae ( at ECL) 0.7356 

Fission troauci poisoning 1.8577 

Burnup considerations 10.2373 

Total inbuilt reactivity .... 30.7863 

i 

j Total Control rod. wort n ..... 34.93 

Burn up (average) - 14657i£ MWa/ton 

anaurance (days of full pswer) - 8^r>f daysn 

! 

Till novi , tne basic reactivity requirements/ 
temperature coeiricitnts a core burnup nave been computed. 
Next tne control system will be evolved ana later certain 
stability analysis none. 
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7.2 KbsC. TC/pCO/i/ob SYSTf.il 

7.2.1 S a lient Features^-; 

Tat.- manoeuvring of warsui^s im r oot.s large 
load Vcxiations on tno reactor. A conventional £v.u control is 
inaat^uate for t;,is purpose as it is aesignea for only a small 
variation from full loau(upto 10%) - tile basic signal it uses 
is tue slow reSj.oriuj.Ly average coolant temperature caanga. In 
marine -reactors/ uowe'/er, an aaoiti. nal prompt signal will be 
required ( Kef (32)) waicn ’feels' tno tovax transient as soon 
as it starts, & aisa^pears graaually at tne ena of loaa variation 
after w.-icn only tile slower reSx-ondiag signal would .oe active. 
Care saoulo be taxn that the reactivity insertion rates ao not 
maKfc tno reactor pronixut critical. For tails / the rate is restricted 
to 0.03 Vo A k (i.o. ax i iroK. 0.04o^) based on the fact that reacto' 
snould respond by a 2.3%/s pawer variation to go from 90% to 20% 
in 30 secs (kef ( 32 ) ) . 

TiiuS/ in the first part of tue operational ; 
tra.usient/ the control system snould act promptly ana efficiently 
performing gross regulation and later fine regulation to prevent 

I 

it ax. x_.ro aci xing unsafe conoitions anu oscillations aro^na 

egu a librium. t 

7.2.2 The Sy^stern : - ' 

I 

'i 

I 

I 

! 

I 

t 

t 
\ 


A simplified uiagram of the system is snown in 



POWER CONTROL 
SYSTEM 



RE ACTOR CONTROL SYSTEM. 










Pig. 7.1. It consists of 2 subsystems, viz. 1 Power control system ‘ 


(prompt acting) and 'temperature control system' (slow acting), 
(i) Foyer Control System:- 


nerc trie main signal is inbalance between 

seconaary steain flow rate, m , and core generated power (proportional 

to neutron flux, 0^) . 2 aaaitioral terms, derivatives of the 

aforesaia, aro aaaed for enhancing prompt cntrol, v/itn scaling 

factors 0. ac fi . Trie rosulciny _ov or error signal activates 
i z 

ei trier control roci cluster vKCC) a, B or C depending on tile signal 
magnitude. " is used to sense traxisient condition; if signal 

exceeds a Ce-rLain veL ue ; D3 (Dv,ad 3ana) is cu&nyed from SP^ 

(wiae, normally) to SP^ (narrow). Rod speed is controlled by an 
error signal similar to tiie above but witn a different gain factor 


b •* The speed control unit assigns a minimum & maximum insertion 
rate ( = 2.3 %/ s j.x)v,er variation) as suuwn in the figure. 


(ii) Ttaupor atur o Control^ 3^ stein s_~ ; 

Here main signal derives from alitor ence between ihe j 

primary coolant average teamr. (T. ;) and a constant reference tempr, j 

av I 

| 

Assent for rod-motion is given wuo-never error signal exceeds tiie | 

• V 1 

tempr. DB. Trie cnar Q ct eristics of the various units are as snown 

| 

in Fig. 7.1. J 

j 

Signals obtained from tne above 2 su.-— systems 
are suaimea up to obtain tna nett error signal. Tne sequential 
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control unit decides tne kcC that is to be moved, with the speed 
assigned by the rod-speed control units. Tne ce.ntrol-rod actuators 
activate tno- particular ROC, txius finally varying 0 ; tnis feedback 
is reflects a in tne power control system resulting in a reduction 
of tne original error signal; tnis effect continues ensuing 
automatic control. 

7.2.3 Control Scheme 

Here, a brief concept of tne roa-coutrol 
strategy over tile core operating lire time is elucidated, 
ti) Control group is constituted of 3 RCC’s A, B anu. C as in 
Fig. 3.3; each cluster moves as a bundle, independent of each 
other. 

(ii) A constitutes the 'Scram' rod, capable of shutting down the 
reactor and is nept fully retracted under normal operation. 

(iii) B constitutes 'coarse' and C, 'fine' control clusters, 

(iv) Initially, B is Jcept fully inserted ana C moved till it 
reaches the top limit. 

(v) Then C is fully inserted and B half extracted; C is moved 
up till it reaches tne top limit. 

(vi) Finally/ C "is fully inserted & B fully extracted; C is moved 
till it reaches the top limit ( Tnis would taen signify tne end 

of core- life). 

Tne above seneme is only a concept & should 
be supported by extensive calculations wnicn is beyond tne scope 
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of tile present work. 


7. 3 _ST^_Ijl. 1TY Are-LYoIS 

7.3.1 Computatiyi'^ of jsom^ parameter sp- 


(a) Tne average neutron life-time, 1* is 


2 2 

1+ L, . 8 

1* = /A 


v 


Eqn (7.15) 


where , v 




a 


= neutron (average) 

= Fuel absorption cross-soc( average) 


1* cnanges witn load as 2, cnanges vjith operating fuel 

3 . 

temperatures. Table 2.3 gives 1* for different loads. 


(b) For neutronics calculations, one-precurs- or group-moael 


nas been cuosen. Tne transier function, G n <s), is given by 


°n < s > = 


r, 


n ( s ) /no 


s+A 


Eqn. (7.1b) 


«r (s) 


1 su (s+A+p/l ) 


wnore, X= 0.013 

p = 0.0065 
Rest of tne 


one-group 
J values • 

symbols have t^eir usual meaning. 


7.3.2 Two - t em pe r at nre Feedback : - 

Tile control loop ror temperature feedback 
effects of fuel and moderator is 
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\(s) 


i V 3 ) 


1(b) 


<X £ .1 £ (a) — 


a__*G- (s) i-s 
m m v ' \ 


6a(s)/n, 



Wt; v<zi< j yi 1 c *z: to ev cili.ia.tfe tne moaer ator . fuel tr an s for functions/ 

G ana G,. 
m f 

Taking an elementary neat balance for a fuel rod in tne reactor 
ana allying a -small perturbation tc tne fuel temper eturo, T^_. 
We arrive at the following relation (from p,92/ ref (34)). 


bgn (7.17) 


l+s'C 


For tne hvoaerator a similar tv^dti. n as aoove is obtoinea:- 


bqn. (7.1b) 


1+ s x m 


where. At = (T f -T c ) at steady -state 

/1\T = (T -T ) at steady state 

m m C 


^ £ = fuel time constant = m^.C^/h^ 
!},-= Power produced per fuel rod 

• A * 

AT f 

t = moaerator time constant. 


Refer appenai^. a 7.1 for evaluation of t ana T , 


l. 


13*3 


Tne 2 f ^ ut'acl jy&ths Ccn oombinuci to onv as 


follows :* 


K ( s ) = 


s+ Y 


Vj'-Ut/ITfe. 


(S-fej-) V2* + ~x~) 

m f 

f. * 

CC -dT 
h = , t f 


hpi (7.1^) 


+ 


a ' At 

m m 


Y = H (c) 


m f 




H ( (_ 1 = c: * At* + a ' ^ 

, h^j f m f + m 


m 


Hence tnc open loop tj_ ansfer function is 


G ( s ) *H ( s ) 
n 


A/1 * (s+ v ) ( s + x ) 


(s+ 1 A iT) ) Cs+ 1/a f ) (o+x-f 3/1 ) 


o-pi (7.20) 


All functions calculated above arc functions of operating 
tamper atur o s and require indiviaual consideration at part- loads. 
Results thus computed are presentee, in Table 7.3. The points of 
intersection of hqn, (7.23) for different loads are evaluated 
by tiie normal rationalizing tocani yie. Calculations are quite 
involved ana uence not snown. The Nyquist Hots consequently 
obtained for aqns. (7.,!^) ana (7.20) are s/etenea in Fig. 7 r 2 
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Table 7.3 Two Temperature Feedback Parameters 


Full 

Load 


£ 


3.693 | 3.409 


3.162 


2.963 


3.693 


. * o 

AT f ( G) 


220.19 156.44 I 99.05 


47.46 


220.19 


(sec) 


0.483 


0.482 


0.481 


0.480 


0.483 


At ( sec ) 
m 


7.72 


5.84 


3.93 


1.98 


7.72 


H (0) 


0.0651 0.0488 0<0326 0.0165 0.0581 


0.0263 0.0206 0.0142 0.0074 0.0255^ 


1.3882 1.4421 1.5088 1.5664 1.2788 


1 (ps ) 11.31 


10.55 


9.82 


9.11 


1-7.43 


Points of Intersection on Imaginary axis of G(s).H(s):- 

w <H ) + 27.74 + 28-02 +28.31 + 28.56 + 22.33 

GH (xl0 ff ) - 6.252 T 6.169 - 6.126 — 6.073 | — • 11.956 

+ + 
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Fie. 7.2 NVQUIST PLOT OF 2-TEMPERATURE FEEDBACK SYSTEM 
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7.3.3 Boiler Transfer Function 



This figure snows 
the lumped model used. 


Assumptions made in 


( i ) Pr imary Temper atur e 


% ( T^+ T, ) 


‘P _ 4 ' bi bo' 

(ii) Secondary Temperature is that of saturated steam, T g 
wnich is a function of pressure. 

(iii) Steam-water mass in secondary side is constant. 

(iv) Throttle is kept constant; power output of boiler is 
proportional to secondary pressure, p . 

using basic heat balance for Boiler Power (BP) the following 
result ensures (from p.99, ref (34)) for Boiler Transfer Function 


® b (s) = K. ■ 4 j ~~ 


Eqn(7.21) 


Where K = 

' x = 


1+ X +S 
X— 1 


X+l 

2 iZ kzZ sJ L. * (1 + 

(T bi- T bo> (T P-V 
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(xn C ) secondary /BP 
t + = -_ B — ■— >- 


V<T bi -V + P„/P 


S' X S 


T = 


(mC ) secondary / BP 
1/< ' r bo- T s > + p/p s 


In the above equation p is the slope of saturation curve , 

s 

(d p /dT ) ' at the appropriate pressure of operation? it is got 

3 

by purturbing the following equation (p in p.s.i. ) 


T (p) = 235.7 + 0.06 12 94 *p - 0.5661E-5*p . 

s 

The secondary side pressure has been chosen as 
500 psi at full ldad & proportionally lesser for part loads thereof 
(by virtue of assumption (IV) above). As before boiler transfer 
functions have to be evolved at various loads and results are 
in Table 7.4* t + and t are still unknown (as secondary side 
design has not been due). The Bode plots, however, for Eqn. (7.21) 

may be plotted, but break frequencies cannot be determined. 

* • _c , , . . , . T , t-- , . which is a design 

As w — * inf, the plots tend to K* T /x + 

parameter^ and at low frequencies, the plots tend to K. 

Fig. 7.3 gives the plots obtained. 

Table 7.4 Boiler Trartfer function parameters . 



Full load 

75% 

50% 

25% 

T 

s 

240.688 

255.732 

204.977 

173.524 

( °C) 

S s 

71.256 

53.255 

48.829 

42.109 

K 

0.7675 

0.9021 

0.9404 

0.9741 

K. T ” 

X + i 

0.628 

0.823 

0.905 

0.964 










138 




ITT 


Fig. 7. 3 Boiler Transfer Function- Body, llot_ 

7 . 3 RESULTS _umD_ CGldiyNTS 

”Tne (Sort Reactivity 3alance is displayed in 

Table 7.2. From figs. 7.2 & 7.3, tne reactor is stable unaer 

all conditions ot operation. Many parameters are However still, 

unknown a re-yiire extensive calculations to be determined- In 

tne reactor control system jz^ and all dead-hands are 

yet to be evaluated & lienee an overall study of tne system would 

not be possible. Transient Analysis & reactor dynamics nave 

not been done in the present study but do require careful 

consideration for c-heck of unallowable overshoots. Throttle 

variations as described in Cnap. 4, ref (34), were studied 

to an extent; but further progress was not possible due to many 

unknown parameters beyond the scope of the present design. 


********************** 



CHAPTER -8 


RESUME 


REACTOR ON SHIP 


8.1 SHIP'S MOTION EFFECTS 

8.2 COLLISION BARRIER & GROUNDING PROTECTION 

8.3 MISCELLANEOUS FEATURES 
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Integration of the reactor with the ship 
includes certLin specific design con siae-reit ions Hitherto not 
discussed lixe snip's motion effects/ induced vibrations, collision 
e grounding protection for reactor, etc. 

o.l Pi'ilj: 1 S woTleis fiCTti On Slsi-ii Ci.iiuasTOjx jt'oRFORH.dsiCh._i — 

Location of tbe reector on board snoula be sued 
as to minimize snip's roll & piten on the steam generator- 
location at about tnu snip’s CG( Centre of Gravity) would be ideal. 

In the- U-tubo steam generator, feea water gains 
neat by a 2-pnase natural circulation system-firstly, feea water 
absorbs tnermal energy from primary coolant becoming partly steam; 
next, tnis steam-water mixture from the tube bundle- passes tnrough 

a steam swirl vane located in the riser whicn separates water I 

i 

i 

from steam cert if ug ally and subsequently through additional 
dr yers to decrease the moisture content. Tnis 2 phase system, 

I 

under the effect. of snips 3-airaen signal motion could present I 

flow oscillations causing loss of heat transfer efficiency, S 

large level fluctuations in the down comer, pulsation in Steam 
Supply, etc., thus forcing a barely stable reactor control 
unit into a region of inst ability . Special simulation computer 
CODfS have to be developed to study the effects of snip motion 
on steam generator circulation; tnis is not attempted here but 
a brief sketch of the associated problem given. 

I 
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The recirculation loop is divided into a number 
of small equal portions of, say, Az. The normal heat balance is 
done for each element. Additionally, the rolling & pitching effects 
subject Az to tangential & centrifugal forces, modifying the 
momentum equation. Tne dipping effect results in a modification of 
'g* (accln. due to gravity) by an oscillatory motion. Consequently, 
the tnermodynamic stability of the recirculation loop in .presence 
of the above effects is done. Finally comes the secondary water 
level control system. Power level, feed water flow & secondary 
pressure dynamically alter tnis level; tnis necessitates the 
development of an auto level control loop to ensure system 
stability. 

Extensive studies in the ship motion effects 
have been done by JM-Sacolo, G.Selvaggi et al (Raf (37)) using 3 
main computer CODES- Sabre- 3M, F-lll MO and DCR-2. The study 
evaluates dynamic behaviour for following perturbations :- 


Roll 

30® (max) 

; period 

= 12 s. 

Pitch - 

10 6 (max) 

; period 

• 

CO 

I! 

Dip 

0. 3g 

; period 

= 7s. 


Results obtained reflect the overall performance of a marine steam 
generator as summed up below:— 

# 

(1) System presents an oscillatory performance for an 
impressed step perturbation. 

(ii) Insertion of an orificed plate suitably positioned in 
the downcomer region increases stability margin & damps the 
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oscillations induced by snip's motion uff^cts. 

(iii) St e am quality increases '..itn< dccrccsin^ water level in the 
downcomer . 

(iv) System stability decreases with increasing steam gtnerator 
inclination or lowor system operating, pressures. 

(v) System stability increases witn lowor pressure crops in heated 
zone or lowor fe_d w,.t c r temperatures. 

B« 2 CGJjLXSIOW ai\IL) GROt: .ib jLikr X ROi'rJCTIe JM 

To minimizo collision off acts for reactor , 
protection is required to avoid reactor aanegu conso. iUo n t i ng in 
uncuntrullca release of fission products (approx. 6 curio-s/W 
generated I ) . T..xs is uone by a grid structure of about SO cm x 
yO cm arranged between longitudnal bulkheads & outerhull, 
supported by decks behind the longitudnal bulkheads (ref (3o) ). 

Tnis structure transfers the collision load into tile transverse 
bulkheads of the reactor compartment, thus destroying the striding 
blow. As per a test program conducted (p.422, ref ibid), no 
mechonical damage will occur for components behind tne collision 
structure; in addition, the protect ien is adequate against 
penetration of a jet aircraft at nominal speed. 

For the greunaing protection, a grid type 
support structure 140 cm high is located above a 200 cm high 
double-bottom (DB). If DB is damaged, the grid is capable of 
withstanding the weight of reactor containment. Consequently, 
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during severe loads , the DB being more easily deformable/ acts 
as an elastic spring. 

8.2.3 MISCELLANEOUS FEATURES : - 

Vibrations are excited by sea-motions/ 
oscillatory propeller forces & auxiliaries; frequency range varies 
from < 1 Hz to 80 Hz. High energy vibrations from propellers are 
specially significant & require detailed analysis. As compared to 
conventional ships of equal power, the machinery weight as well as i 
the normal cruising load on machinery is higher by 1.7 & 3.5 times j 
respectively (approximate figures). Tne Power consumed by the ! 

\ 

primary coolant pump is substantial. ' 

The ship 1 s speed depends on the SHP ( Shaft 
Horse Power) developed by the engine. The speed versus SHP 
characteristics of a ship is dependent on the displacement, profile 
at the water line & certain other factors; it however follows an 
approximate cubic law with linear speed increments requiring cubic 
SHP increments. Consequently, at 50% of power developed, the 
ship acnieves nearly 75% of its maximum speed* A rough estimate of* 
22 knots may be expected as the maximum speed of the ship presently 
under consideration, with a 30,000 SHP maximum power output ( assuming 
70% of steam to be used for propulsion only, and a plant efficiency 
of 22%). 

An important aspect hitherto ignored is 
requirement of a secondary propulsion system. Coupled to the main 
gear box, a compact gas turbine of approximately 15000 SHP may be 
included so that, in the event of a main engine breakdown, the ship 
is still capable of making way at 65% of its maximum speed which, 
now, would correspond to that of 45000 SHP (**25 knots). 



CHAPTER -9 


CONCLUSION 


RESUME 


9.1 SUMMARY 


9.2 SUGGESTIONS FOR FURTHER STUDY 
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By now, the salient aspects of marine 
reactor design have been covered. In conclusion, the destroyer 
class ship is about the smallest to which nuclear propulsion can be 
reasonably applied - a preliminary study of containment weight revea 
this (chap 6), as a propulsion machinery weight of ^ 2500 tons is 
bound to accrue ( & a displacement of the ship of ^8000 tons). 

9.1 SUMMARY 

Some of the design parameters evolved are ' 

t 

briefly summed up in Table 9.1. Sinc^ tne various design processes 

! 

are inherently interactive, change of a parameter may have a ' 

retrospective effect in other processes. In the shielding design [ 

I 

a machinery space of 300 cms (instead of the 200cms assumed) would j 
have resulted in a dose rate ofo.4982 rnrom/hr outside the 
containment with a weight increase of 262.0 tons and diameter 
of 6.6 ft. This space can be fixed only after all components of 
primary/secondary loop are designed. Core burnup obtained has been 
low for a 5% enrichment. Tnistimmalation is due to the 

compact, close— latticed core configuration cnosen that is so 
vital in reducing overall containment weight in marine reactors. 

Table 9.1 Design Characteristics 

I GENKRAD REACTOR DATA 

A. Reactor Thermal Power .... 120 MW 

No. of Primary coolant loops ... 2 
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3. Primary Coolant 

- System pressure ...... 2000 psi 

- Total flow rate 5553 tons/hr 

- Inlet tempr. 280 °C 

- average rise in core 15.42 °C 

C. Fuel (V0 2 ) 

- No. of core regions — 1 

- No. of fuel asslys. — - — 29 

- No. of fuel rods per assly. — 282 

- Enricnment 5 w % 

— Core Power density .... 63.51 MW/m 

- Eqbm. core average burnup ... 14,697.8 MWd/t 

- Wt. of U0 2 in core .... 6.919 tons 

D. Control rods 

m 

— Type/Number RCC/29 

- Absorbing Material B^C 

- Rods outer dia/pellet dia ... Qr 51/0 ,494 cm 


E. Secondary System 
- Steam pressure (at full load) 

•f Steam quality .... 

F. Plant Physical Parameters 

— Overall Wt. of shielding 
— Containment outer, dia 

II. NUCi£AR DESIGN DATA 

A. Structural characteristics 

- Core dia/height .... 134.29 cm/133. 4cm 

- Fuel rod & control rod cladding .... ss347 

- Fuel cell type/pitch ... square/1. 30cm 

- Fuel rod dia/pellet dia ... 0.99/0.858 cm 


500 psi 
0.9975 


838.2 tons 
30.67 ft. 
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-Volume fraction H 2 g/U 0 2 in core ..... 1.592 

B. Pnysical & performance cnaracteristics 

- Fuel type/density sintered U0 2 /10. 97gm/cc 

- Fuel Enrichment - 5 w% 

- Hot channel factors , 

F A J / F% /total F. 3.44/1*04/3.58 

U Q w 

F*J / F e u / total F „ 1.994/1.22/2.105 

R H ri 

- Flux functions analysed aSin(z)/Sin(z)/(w-z)* 

BOL/Normal/EOL Sin(it-z) j 

C. Control characteristics 

_ Type mechanical/ RCC | 

- Absorbing material .... B^C(den 0.7 gm/cc) 

- Rods per assembly /No .of . asslys. ... 6/29 

D. Kinetic characteristics 

— Moderator tempr. coeff. (full load) ...-6.233E-4 

Ak/k 

-Fuel tempr. coeff. (full load) . . . . -2. 738E-4 £k/k 

III. Core Reactivity Balance 
Refer Table 7. 2 . 

9.2 Suggestions for further study 

Some important aspects of marine reactor 
design not done in the present work are : - 

(a) Materials design -■ The performance of only a few materials 

has been studied; optimization study of different reactor materials 
is suggested to evolve better design parameters. 

(b) Ship-motion effects on the steam generator require a study 
in far greater depth than mentioned in chapter 8. This proolem 
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is peculiar of marine reactors and has a direct effect on overall 
system stability. 

(c) In-core instrumentation design requires careful consideration 
as it conveys core fluxes and temperatures at various points in 
the core; these signals structure the basic inputs for the reactor 
auto-control loop. Control rod wor'th evaluation is also suggested. 

(d) In the Snieldiny Design the radiation from reactor effluents 
requires evaluation after steam-generator and piping design is 
done; its magnitude is, however, small & tne shielding thicknesses , 

I 

may require marginal modification to counteract this additional dosf 

(e) The pressure-vessel & piping design for mechanical, thermal 

and other stresses necessitates careful checking out; in the | 

present design a vessel thickness of 12 cm has been assumed I 

| 

based on available literature and not on calculations. | 

(f) Design of pressurizer, steam generator, components of j 

t 

secondary loop & subsequent heat balance is also proposed. 

(g) Reactor start-up & shutdown procedures are suggested for 
the present design. 

Study of an extensive topic as a reactor 
design in a limited time frame effects in simplifications & 
exclusions. The core flux functions for thermal hydraulic study 
do not precisely represent actual flux variations. Other flux 
functions as considered adequate may be assumed to evolve more 
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realistic performances. The core nuclear design (Reactor Physics 
& Shielding ) has been done using a heterogeneous 2-group theory 
with constant core temperatures for fuel & moderator- 
multigroup theory using local temperature variations and 
development of transport theory computer CODES is suggested; 
this evolves a more precise design , closer to a practical 
reactor unit. The reactor automatic control loop may be studied 
in detail, both for system stability & transient characteristics, 
once the loop constants (of Fig. 7.1) are evaluated. Study 
of reactor transient for an ahead-astern simulation consisting 
of a load variation from 90% to 20% in 30 secs and immediate 
return to 90% in another 30 secs would be of significance in 
ascertaining the system performance. 

An earnest attempt has been made to describe 
the complex, varied & multi-dimensional problem that a warship 
reactor design is. Both strategically and tactically, nuclear 
propulsion is the ultimate form of propulsion for warships 
today; it is obvious to all who pause to ponder this matter. 

It will remain so until a superior form of ocean propulsion 
for supporting modern weapon systems is architectured. 
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— (zSinz) =0 at maximum ; 
da 

•=> Tan(z) = -z 

The above eqn. is satisfied when z=2. 02877 
So, maximum = 2.02877*Sin- (2.02877) = 1. 8197 


We shall choose the axis as above in diagram. 
H 


x av 


H 


r z. sj 

w 1.8 


Sin.z.dz 


where, e= 


H -H 
e 


8197 


As z is distance in cms we will have to convert it to radians* 


= % 


H 

rv 


q av - -y * 

H ^ 


^ "t 

* Sin( — ) *dz * — 1 


H 


e e 
% _ 


H 


^dx 


TC 


H 


.dz 


1.8197 


Let 


x 



15b 


1.8197 * H 


x. Sin x.dx 


= 2o 

1.8197 


I. 

v H ' * % 


; hh _ ( ,h % ) + 


Cos (- 


-) + Sind 


Sin (- 


Putting A = 


B _ eJL 

' B ~ H 


(H -H) 71 
e 


N 

F„ - **0 


1.8197* ^ * H 

4 ... 

H e (-ACos(A) + BCos(B) + Sin(A)-Sin(B) > 
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3. TRUNCATED BESSEL* S FUNCTION : 



Function 


F(r)=AJ Q (ly) + CI Q (mr) 
R = core radius 


A,C, 1 & m are all constants determined 

from trie Reactor Physics CODE and hence known. By definition 
H R 

1 


av -m 2 H 


J7 

o o 


F(r) * (Axial functipn) *2 rdrdz 


The axial function integration is done in (1) and 

( 2 ) above . 

The radial integration is:- 
R 

q &v = f * A- J Q (lr ) + C- Jq (mr ) ) * 2.r.dr 

R 0 

Since j Jlr . Jq (lr ) .dr = lr . (lr) 


and 


s» 

mr*J^(mr)dr = mr * 1^ (mr)/ 


%V =^2 ( f--R-a i( iR) + £-. R .x i(nR) , 
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At r= 0 , the peak flux occurs; its value is (A+C) 


^av 


A + C 


2'A- J^U-%) + 2 C • I 1 (mR) 


1R 


mR 


1 + 


A 


2 . J <1R) 


1R 


+ 2( ~) , I 1 (mR) 


mR 


As 


1 / 


N 


2 J (1R) 


1R 


+ 2 - ( f- )• J l (mR) 


mR 
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A 3.2 DETERMINING K : 

Let K fu <T) - a+ bT ^ be a linear function of temperature. 

AIM To determine what value of 

K- (T) is to be usect in the 
fu 

heat conduction equation, 

Eqn. (4.3) . 


Solution 

We shall develop an equation relating T Q and T g 
using K (T) ~ a+bT then compare tnis with Eqn»(4.3). 

In the steady-state condition 


.%-r 2 Al = - K fu (T).2'jrr.A 1. , in any fuel ring 

of thickness dr at radius r. Integrating from r=0 to r— r fu 


. ' ' ' , r 2 , r fu 

* q < o~ ) 

2 0 

i 1 1 

. . q 

(or) - 




*■ - 2 <aT *§£ )l £ 
z 0 

r 2 


T* 4-T 

4 ( a+b : ( 0f j) 


This means (T) in Eqn (4.3) is evaluated at the arithmetic 

T + T 

.mean fuel temperature for determining i.e. at T = _0 f_ 


2 
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A 5.1 CALCULATION OF D &g i- 


Element 

NUM DEN* 

( 

2a* 

a s' ' 
0>) 

: 

( 1- P-q) 

2* 

F 

e 

6.044 E22 

0.093180 

11.4 

0.68902 i 

3. 9880 

0.68075 

C 

r 

1.6719E22 

0.030758 

4.2 

0.07022 i 

3. 9871 

0.06931 

N. 

1 

8. 2291E21 

0.025058 

17.0 

0.139895 

0.9885 

0.13829 

s. 

1.7536 E21 

0.000178 

2.2 

0.003858 

0.976 

: 

0.00376 

N b 

4.1611 E20 

0.000290 

6.5 

0.002705 

; 0.9928 

0.00269 

M 

n 

1.7541 E20 

0.001401 

1.9 

0.000333 

0.9878 

0.00033 

C 

i 

3.2091 E20 

j 0.0000009 

4.7 

0.001500 

0.944U 

0.00142 

| 

Total 

0.150867 

! 

0.90754 


3.896555 


* Taken from Table A6 . 1 


3(£ a +£ tr> 


0.318241 cm 
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A 5.2 CALCULATION OF D 


The fast diffusion constant for core, is 


, given by the expression :■ 


Eqn. (1) 


lf ^<2 la+ ^a-fr 0 )) 

where, ^ia = av “ fast absorption; cross section. 

= averaged scattering cross section, 
since, the scattering cross section of Hydrogen varies 
strongly with neutron energy, we use an artifice to calculate 
its scattering cross-section. Let us make an initial guess 

of = 10 barns - 

The parameters are 



H. 

0 (total) 

S. S. 

Total 

11 I 

n ( #/cm 3 ) 

0.3641E23 

0.3622E23 

0. 8407E22 j 


% (b) 

10 

CM 

• 

11.4 



1.0 

0.$2 

0.353 


U-iy 

0i 333 

0.258 

0.988 


7 , 

s 

0.5094 

| 

0.15212 

0.09584 

0.7577 


For the above parameters. 




(tot) 


a 


0.6862 
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(1 


-pj 


i-i 


TXjB.C 


2 


r = 

s 


ln( §». ) 


S (TOT) 

ft 


= 0i.5520S 

h ' 


1 

* tr 


0.7475 


3 r |-<i- h > 

We see above & one from tabulation do not agree. 

So we guess hicher values of & s (H 2 ) and continue calculations 
We obtain convergence when <§T (h ) = 13.7 

3 2 

"^roac ~ 0.69693 


(1 “^nac 

2 


= 0.54459 


= 0.7488 cm ” 1 Eqn. (2) 


aioac 

Next we have to determine 2 la . 


•#1 


N" I 


efP 


ln( _ 0 _ ) 
E th 


we 


know from calculation for 'p 1 , that 


1&f£ A+ ^S 


= 21.254 barns for U0, 


2 ia = 2.6323E22 * 21-254E-24 = 0.03858 cm 


-1 

Eqn (3) 

From Eqns. ( 2 ) & (3), we can evolve D-^ "through Eqn (l) 


Dt, 

rx 


3 * (0.03858 + 0.54459 * 0.7488) 
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AS. 2 INCOMPLETE ERROR FUNCTION* 


JLt I JL 


40 


n 


j 


- tx 


-xn 


it 


<3t = x 


n-1 


-t 


x 


Lll 


dt 


are called incomplete error functions. 
They nave series solutions. 


E Q (x) = e“ X /x 


E n (x) - - E n-l (x) 


In general/ 


<-x ) 


m 


m*o 

nufn-1 


(n-l-m)m .* 


n n-1 

+ (-1) 7^.1) j<V +ln(x)-A n ) 


S n (x) “ X 

( n— . f —m im * — — . - 

^ n>- o 
n-1 

where/ f = 0.577216 / A ± = 0 , A n = ^ ( “) 

■*n=l 


In particular 

E x (x) = - ( y +ln x) + x- x 1 x 

2 3 

E 2 (x) = 1+ x(f+ In x) - + jj 


32 


The following integral has been used i- 

n- 

X 


J kx n-1 i kx 

e .E n (ax+b).dx = ^ ~^J5T~’ e % E - 4 (ax+b) - 


n-i' 


i=0 


n-1 bk . a* 4 ** . 

a a”. E, ( p — ) 

“7n“*e * 1 1~ *- 

k a 


Ref ( 31 page 355 
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In particular. 


Joe 

e . ( ax+b ) . dx 


> 10 . 


“ — E (ax+b) - e bk//a * _ ■ ax+b , 
k 1 k E l ( 1-k/a > 

E r (x) is well approximated by:- 


E (x) = e“ x Q (x) 
n n 

1 n 

where, Q (x) =— — + — - — r 

n x+n (x+n; 


A6.3 


(a) Solution ; ^ 0.(r) - c 2 *0. (r) + _ Qe~^ = 


-2r 

V 2 0 s (r) - c 2 * 0 s (r) + P % 


In spherical coordinates 


= -i— * ( x: 2 ) 

v 2 dr ^ r * dr ; 


where P = Qe 


r < + 2?_fs_ > - = 2 * <*> + 

dr dr s r 


Substituting r _(r ) - X’ • equation reduces to 


<D 2 -c 2 ) X = - p £ 

r 


- r2!v 


By using the simple integrating factor method by first 
considering (D+c) X=Y , solving for Y and later, 
solving for X, we get the equation 
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-2cr 

X = + Be cr 

-2c 

From appendix A6.2, 
As x = 0 g *r 


cr f -2cr f e <C ~^ )r 
- P.e J dr.c ^ cr J ~ — .dr 

. cl" - - E 1 < (c-2 ).r ) 


/ * r 


Ae 


-2cr 


Be 


cr 


2c 


+ Pe 


“Jdr. 


-2cr 


■B^r, (<*£)) 


Again using the integral identity of A6.2, 


^ s (r) - - Hr 


-2cr 


Be 


cr 




cr 


(c-S>) 

2c 


r l 2c ' E 1 ( “ 14 2 c_ 

v C— d 

-^iP” *V (c - 2)r > < 


(c-Z 


c-r 


AS 


►- 0 as r — vinf , B = 0 is essential. 
Thus, substituting back for P, we obtain:- 


_2cr Xj < y >7 . / \ 

0 s ( r) = _A|__ + ^( c -2)*E i (Y)e Cr “ E l(X)*e _Cr j 

verier e, X = (c-Z).r 

Y = >7(1+ ^ ) 

(b> Solution of 0 (r) - c 2 . 0 (r) + Qe~ Er = 0 
— . s s 

Here the approach is similar. Substituting for 


r * 0^(r) = X, we obtain :- 
<D 2 -c 2 ) X = -Pr.e“ r ^-. 


v 


Eqn (1) 
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The complementary Function of -i iis differential Eqn. is 

-cr cr 

X = Ae + Be 

The particular Integral is obtaineo by choosing a particular 
solution for X = (Dr + E)e~ ; 


Putting this is eqn. (l), we obtain the coefficients 
1 . ™ _ 2 


D = 




i E 


( c 2 - Z 2 ) 2 ' 


-Zr 


-*CJT cr 

Hence , *5_ ( r) = + Q e ( r - 2 ) 

3 r r r cTZp (c2 _£)*- 

as & (inf) =9, B=0 
s 


• # (r) = ^ 

s r 


-cr 




2 .i 


(c 2 -! 2 ) 


_ ) 
2 ) 


★ **★**★* *‘* ★★**■*****•*■*’* ,v 


TOTAL GAMMA EMISSIONS FROM COKE 
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A 7.1 EVALUATION OF MODERATOR TIME CONSTANT {*& ) 


Heat balance for a fuel rod channel gives 
teat added) = (heat suppliec to co< Lant) - ( heat carried away- 
flowing coolant ) . 


r) mCp d T ac 
P dt 


Ah (T ,~T ) + ih L (T.-T ) .... Eqn(l) 

c cl ac p x ac 


where, Tj_ , T = coolant inlet/average teinprs. in core 
ac ~ ^ 

T ^ = c ladding outer t empr. 

A = coolant flow rate 


jninatiny T cl in Eqn(l) by substituting 

got by Steaoy-state 


h. T , + h T 
t - £ af c ac 

cl h^ + h 


t transfer between fuel & coolant, where h^, T ^ are fuel 
- transfer coefficient and average tempr. given byEqn(7. 17) / 
>btain:- 


d T . 
ac 

dt 


A h 


^e^f ^ mC p 


( T — T ) - — 

af ac in 


■(T -T. ) . . Eqn<2) 

ac x 


e moderator time const-ant is the coefficient of its differential 


m ~ <Jfa*VGa c *h f >aCh j” 1 - £qn(3) 

; evaluated from Eqh(2) by' putting dT ac 

dt ~ ° 

tr Steady State condition). We obtain 


1 



Eqn (4) 


h 

c 


1 4. 

_ — V 1+4DH £ 
D 



where, D = A- (T 


a£ 


) / rn C 
ac' / ^ 


(T -T. ) 
ac x 


since all terms in Eqn ( 3 ) are 
temperatures change with load 
load variations as -well 


kncv,*n,T;^ is evaluated. As 
variations, varies with 



